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Abstract 

Given n independent random marked d-vectors Xi with a common den- 
sity, define the measure = X^j^i, where is a measure (not necessarily a 
point measure) determined by the (suitably rescaled) set of points near Xi. 
Technically, this means here that stabilizes with a suitable power-law decay 
of the tail of the radius of stabilization. For bounded test functions / on R'^, 
we give a law of large numbers and central limit theorem for I'nif)- The latter 
implies weak convergence of fn(-)) suitably scaled and centred, to a Gaussian 
field acting on bounded test functions. The general result is illustrated with 
applications including the volume and surface measure of germ-grain models 
with unbounded grain sizes. 
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1 Introduction 



This paper is concerned with the study of the hmiting behaviour of random mea- 
sures based on marked Poisson or binomial point processes in d-dimensional space, 
arising as the sum of contributions from each point of the point process. Many 
random spatial measures can be described in these terms, and general limit theo- 
rems, including laws of large numbers, central limit theorems, and large deviation 
principles, are known for the total measure of such measures, based on a notion of 
stabilization (local dependence); see [16, 17, 18, 21]. 

Recently, attention has turned to the asymptotic behaviour of the measure itself 
(rather than only its total measure), notably in [2, 3, 9, 14, 19, 21]. It is of interest 
to determine when one can show weak convergence of this measure to a Gaussian 
random field. As in Heinrich and Molchanov [9], Penrose [14], one can consider a 
limiting regime where a homogeneous Poisson process is sampled over an expanding 
window. In an alternative limiting regime, the intensity of the point process becomes 
large and the point process is locally scaled to keep the average density of points 
bounded; the latter approach allows for point processes with non-constant densities 
and is the one adopted here. 

A random measure is said to be exponentially stabilizing when the contribution 
of an inserted point is determined by the configuration of (marked) Poisson points 
within a finite (though in general random) distance, known as a radius of stabi- 
lization, having a uniformly exponentially decaying tail after scaling of space. This 
concept was introduced by Baryshnikov and Yukich [3], who proved general results 
on weak convergence to a limiting Gaussian field for exponentially stabilizing mea- 
sures. A variety of random measures are exponentially stabilizing, including those 
concerned with nearest neighbour graph, Voronoi and Delaunay graph, germ-grain 
models with bounded grains, and packing; see [3]. 

In the present work we extend the results of [3] in several directions. Specifically, 
in [3] attention is restricted to the case where the random measure is concentrated at 
the points of the underlying point process, and to continuous test functions; we relax 
both of these restrictions, and so are able to include indicator functions of Borel sets 
as test functions. Moreover, we relax the condition of exponential stabilization to 
power-law stabilization. 

Along with the central limit theorems, we give a law of large numbers for random 
measures (under weaker stabilization conditions) which is proved using similar ideas. 
In Section 7, we illustrate our general results with applications to germ-grain models 
with unbounded grains, and to questions of Voronoi coverage raised by Khmaladze 
and Toronjadze [10]. Many other fields of application have been discussed elsewhere 
[2, 3, 16, 18] and we do not attempt to review all of these. 

We state our general results in Section 2. Our approach to proof may be summa- 
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rized as follows. In the case where the underlying point process is Poisson, we obtain 
the covariancc structure of our limiting random field (and also to prove the law of 
large numbers) by a refinement of the approach to second moment calculations used 
in [3] and [18], often known as the objective method [1, 22], which is discussed in 
Section 3. To show that the limiting random field is Gaussian, we borrow normal 
approximation results from [19] which were proved there using Stein's method (in 
contrast, [3] uses the method of moments). Finally, to de-Poissonize the central 
limit theorems (i.e., to extend them to binomial point processes with a non-random 
number of points) , in Section 6 we perform further second moment calculations us- 
ing a version of the objective method. This approach entails an annoyingly large 
number of similar calculations (see Lemmas 3.8 and 6.1) but avoids the necessity 
of introducing a notion of 'external stabilization' (see Section 2) which was used to 
deal with the second moment calculations for de-Poissonization in [3]. This, in turn, 
seems to be necessary to include germ-grain models with unbounded grains. We 
give our proofs in the general setting of marked point process, which is the context 
for many of the applications. 

2 Notation and results 

Let {Ai, J-'m, fJ'M) be a probability space (the mark space). Let d E N. Let 
^(x, t; X, A) be a Borel-measurable M- valued function defined for all 4-tuples (x, t, X, A), 
where X cW^ x M. is finite and where (x, t) e A" (so x G M°' and t G M), and A is 
a Borel set in R*^. We assume that ^(x, t; X) := ^(x, t; X, •) is a measure on R'' with 
finite total measure. (Our results actually hold when ^(x, t; X) is a signed measure 
with finite total variation; see the remarks at the end of this section.) 

Suppose (x, t) G M"' X TVf and X dW^ x Mis finite. If (x, t) ^ X,we abbreviate 
notation and write ^(x, t; X) instead of ^(x, t]X U {(x, t)}). We also write X^'^ for 
A'U{(x,t)}. 

Given a > and y G M*^, we let y -|- aX :— {(y -|- ax, s) : (x, s) G X}; in 
other words, scalar multiplication and translation act on only the first component 
of elements of x A4. We say ^ is translation invariant if 

^{^,t;X,A)^^{y + K,t;y + X,y + A) 

for all y G M*^, all finite X C M."^ x M and e X , and all Borel A C M'^. Some of the 
general concepts defined in the sequel can be expressed more transparently when ^ 
is translation invariant. 

Let K be a probability density function on W^. Abusing notation slightly, we 
also let K denote the corresponding probability measure on W^, i.e. we write k,{A) 
for «;(x)(ix, for Borel A C M"^. We shall assume throughout that the density 
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function k is bounded with supremum denoted ||k||oo, and that k is Lebesgue-almost 
everywhere continuous. Let supp(«;) denote the support of i.e., the smallest closed 
set B in with k{B) = 1. 

For all A > 0, let V\ denote a Poisson point process in M.'^ x A4 with intensity 
measure Xk x /im- For n let Xn be the point process consisting of n independent 
identically distributed random elements of x Al with common distribution given 
by K X fij^. Let Hx denote a Poisson point process in M*^ x Al with intensity A times 
the product of d-dimensional Lebesgue measure and ^j^. 

Suppose we are given a family of Borel subsets Fa of M'^, indexed by A > 1. 
Assume the sets F are nondecreasing in A, i.e. Fa C Tx' for A < A'. Denote by F 
the limiting set, i.e. set F = Ua>iFa. We assume that the topological boundary 
of F (i.e., the intersection of the closure of F with that of its complement) has zero 
Lebesgue measure, and that int(F) = UA>iint(FA), where int(-) denotes interior. 
The simplest special case has Fa = M'^ for all A. 

For A > 0, and for finite X cW^ x M with (x, t) e X, and Borel A C W^, let 

eA(x, t; X, A) ^(x, x + AV'^(-x + A"), x + + ^))lr,(x). 

When ^ is translation invariant, the rescaled measure ^a simplifies to 

eA(x,t;A',^) =e(A^/'x,t;AV'^A',AV'^A)lr,(x). (2.1) 

Our principal objects of interest are the random measures /x^ and z/^^^ on IR*^, 
defined for A > and n e N, by 

4:= ^A(x,t;PA); 4„:= ^A(x,t;A'„). 

(x,t)GPA (x,t)GA'„ 

We arc also interested in the centred versions of these measures JZ^ := //^ — E [/i^] 
and T^xn '■~ ~ ^[^\n\ (which are signed measures). We study these measures 
via their action on test functions in the space BiW^) of bounded Borel- measurable 
functions on M"*. We let B{M.'^) denote the subclass of B{M.'^) consisting of those 
functions that are Lebesgue-almost everywhere continous. 

Given / e S(R'^), set {f,U^,t-X)) := j^, f {7)^^,1- X ^dz). Also, set 




and set (/,7li) := J^cfdfii, so that {f,J]\) = (Z,/^!) -E{f,^\). Similarly, let 
(L 4,n) fd4,n and := (/, - E (/, u^. 
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Let I • I denote the Euclidean norm on W^, and for x e R'' and r > 0, define 
the baU Br{x.) := {y G M*^ : |y - x| < r}. We denote by the origin of IR*^ and 
abbreviate -Br(O) to Br- 

We say a set A" C M"^ x is locally finite ii Xn{B xM) is finite for all bounded 
B C W^. For (x, t) e X M and Borel A C M'^, we extend the definition of 
^{x,t;X,A) to locally finite infinite point sets X by setting 

^{x,t;X,A) := limsup^(x,t; [Xn{BK x M)],A). 

Also, we define the x-shifted version {^*(-, ■) of ^(x, t; ■, ■) by 

e^*(A',A)=e(x,t;x + A',x + A). 

Note that if ( is translation-invariant then ^^\X,A) = ^^\X,A) for all x e R'^, 
teM, and Borel A C W^. 

Notions of stabilization, introduced in [16, 18, 3], play a central role in all that 
follows. 

Definition 2.1 For any locally finite X gW^ x M, any (x, t) e R'^ x M, and any 
Borel region A C R'', define R{yi,t] X , A) (the radius of stabilization of ^ at (x, 
with respect to X and A) to he the smallest integer-valued r such that r > and 

e(x, t;^+{[Xn{BrX M)] uy),B)= e(x, t;^+[Xn {Br X M)], B) 

for all finite 3^ C (A \ B^) x M. and all Borel B QW^. If no such r exists, we set 
R{x,t;X,A) = oo. 

When A is the whole ofR'^ we abbreviate the notation -R(x, t; X, R'^) to i?(x, t; X). 

In the case where ^ is translation-invariant, i?(x, = R{0,t;X) so that 

i?(x, t; X) does not depend on x. Of particular importance to us will be radii 
of stabilization with respect to the homogeneous Poisson processes Tlx and with 
respect to the non-homogeneous Poisson process Vx, suitably scaled. 

We assert that R{^,t] X , A) is a measurable function of X, and hence, when X 
is a random point set such as Tlx or Vx, R{^, t', A) is an N U {cxD}-valued random 
variable. To see the assertion, observe that by Dynkin's pi-lambda lemma, for any 
A; e N the event {R{:s.,t] X , A) < k} equals the event nBeB{s{X,B) = i{X,B)}, 
where B is the Il-system consisting of the rectihnear hypercubes in R*^ whose corners 
have rational coordinates, and for S e B we set 

s{X, B) := sup{e(x, t-^+ {[X n{BkX M)]Vjy),B) -.y (Z {A\ B^) x M], 
i{X, B) inf{^(x, t; x + {[X n {Bu x M)]uy),B) -.y Q {A\ Bk) x M). 
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Also, s{X , B) is a, measurable function of X because we assume ^ is Borel-measurable 
and for any b we have 

{X : s(X) >b} = 7ri({(A', y) : ^(x, t;^+[Xn{Bk x M)] U [y\{B, x M)], B) > b}) 

where tti is projection onto the first component, acting on pairs {X, y) with X and 
y finite sets in x A4. Similarly, i{X, B) is a measurable function of X. 

Definition 2.2 Let T, T' , T" and T'" denote generic random elements of M. with 

distribution fiM> independent of each other and of all other random objects we con- 
sider. Similarly, let X and X' denote generic random d-vectors with distribution k, 
independent of each other and of all other random objects we consider. 

For p > 0, we consider ^ satisfying the moments conditions 

sup ¥.[U-^,T-VxXY]<^- (2.2) 

A>1, xesupp(/t) 

and 

sup E[eA(x,r;P;,U{(y,T')},M"r] <oo. (2.3) 

A>1, x,yesupp(«;) 

Let ^*(x, t;A',-) be the point measure at x with the same total measure as 
^(x, t- X, •), i.e., for Borel AQW^ let 

r(x, t- X, A) := e(x, t- X, M'^)Ia(x). (2.4) 

We consider measures ^ and test functions f & B (M*^) satisfying one of the following 
assumptions: 

Al: ^ = i.e., ^(x, t; X, ■) is a point mass at x for all (x, t, X). 

A2: ^(x, t] X, ■) is absolutely continuous with respect to Lebesgue measure on R*^, 
with Radon-Nikodym derivative denoted ^'(x, t;A',y) for y G K*^, satisfying 
^'(x, t] X, y) < Kq for all (x, t; X, y), where Kq is a finite positive constant. 

A3: / is almost everywhere continuous, i.e., / e S(R'^). 

Note that Assumption A2 will hold if ^(x, t, A",-) is Lebesgue measure on some 
random subset of M*^ determined by x, t, A'. 

Our first result is a law of large numbers for (/, /ij) and (/, i/^ ^), for / e S(M'^). 
This extends a result in [18] where only the case where / is a constant is considered. 
We require almost surely finite radii of stabilization with respect to homogenous 
Poisson processes, along with moments conditions. 
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Theorem 2.1 Suppose that i?(x, T; ^^(x)) is almost surely finite for K-almost all 
X G r. Suppose also that f G B{R.^), and one or more of assumptions Al, A2, A3 
holds. Let g = 1 or q = 2. Then: 

(i) If there exists p > q such that (2.2) and (2.3) hold, then as X ^ oo we have 

A-'(/,/xi) ^ ^/(x)EeS;^(7i,(x),M'^)«:(x)rfx (2.5) 

and 

J2 \{f,U^,t;Vx)-Cx{^,t;Vx))\^0. (2.6) 
(x,t)ePA 

(a) If \{n)/n 1 as n —>■ oo, and there exists no > and p > q such that 

sup E [Cm„){X, T; Xn-iXT] < oo, (2.7) 

n>no 

then we have as n ^ oo that 

r^~\fAin),n) ^ jj{^m^o^\n.i.)X)<^)d^. (2.8) 

and 

Yl \{f,U^,t;Xn)-Cx{^,t;Xr,))\^0. (2.9) 

Our main results are central limit theorems to go with the laws of large numbers. 
For these we need further conditions. The first of these requires finite radii of 
stabilization with respect to homogeneous Poisson processes, possibly with a point 
inserted, and, in the non translation-invariant case, also requires local tightness of 
these radii. 

Definition 2.3 For x G M'^ and \ > 0, we shall say that ^ is A-homogeneously 
stabilizing at x if for all z G M*^, 



lim sup max{R{y,T;7i\),R{y,T;T-C^ )) < oo 



1. (2.10) 



In the case where ^ is translation- invariant, i?(x, t;A') does not depend on x, and 
^^^(•) does not depend on x so that the simpler-looking condition 

P[R{0, T; Hx) < oo] = P[R{0, T; Tff) < oo] = 1, (2.11) 
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suffices to guarantee condition (2.10). 

We now introduce notions of exponential and power-law stabilization. The ter- 
minology refers to the tails of the distributions of radii of stabilization with respect 
to the non-homogeneous point processes V\ and Xn. 

For A; = 2 or A; = 3, let Sk denote the set of all finite A C supp(fi;) x M. with 
at most k elements (including the empty set), and for nonempty A & Sk, let A* 
denote the subset of supp(K) x M. (also with k elements) obtained by equipping 
each element of A with a /i^^-distributed mark; for example, for A = {x, y} e S2 
set ^* = {(x,r'),(y,T")}. 

Definition 2.4 For x G W^, A > and G N, and A G S2, define the [0, oo\-valued 
random variables Rx{x,T) and Rx^n{^,T; A) by 

i?A(x, T) = i?(x, T; AV'^(-x + Vx), A^/'^(-x + supp(«;))), (2.12) 
i?A,n(x, T; A) = i?(x, T; V/'^(-x + {X^ U A*)), A^/'^(-x + supp(«))). (2.13) 

When A is the empty set we write i?A,n(x, t) for -R(x, t; 0). 

For s > anc? e G (0, 1) define the tail probabilities t{s) and Ts{s)by 

t{s)= sup P[Rx{^,T) > s]; 

A>1, xGFa 

T.(S)= sup P[i?;,,„(x,r;^*)>s]. 

A>l,nGNn({l-e)A,{l+e)A),xer;^,y4e52 

Given q > 0, tt;e say ^ is power-law stabilizing of order q for k if sup s'^r(s) < 00. 
We say ^ is exponentially stabilizing for n i/ limsup^^^o log r(s) < 0. We say 
^ is binomially power-law stabilizing of order q for k if there exists £ > such that 
sup5>i s^T^{s) < 00. We say ^ is binomially exponentially stabilizing for k if there 
exists £ > such that limsup^^^ logTe(s) < 0. 

It is easy to see that if ^ is exponentially stabilizing for k then it is power-law 
stabilizing of all orders for k. Similarly, if ^ is binomially exponentially stabilizing 
for K then it is binomially power-law stabilizing of all orders for k. 

In the non translation-invariant case, we shall also requre the following continuity 
condition. 

Definition 2.5 For x G W^, we say ^ has almost everywhere continuous total 
measure if there exists Ki > such that for all m & N and Lebesgue- almost all 
(x, xi, . . . , Xm) G (M'^)™'+^, and {hm x ■ ■ ■ x iJiM)-cdfnost all {t, ti, ^2, • • • , tm) £ M"^'^^, 
for A = M*^ or A — Bk with K > Ki, the function 

(y, yi, y2, • • • , ym) ^ C(y, t; {(yi, ti), (y2, ^2), • • • , (ym, tm)}, y + A) 

is continuous at (y, yi, . . . , y„) = (x, Xi, . . . , x^). 
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Subsequent results will require one of the foUwing formal assumptions to hold. 
A4: ^ is translation- invariant. 

A5: ^ has almost everywhere continuous total meaure. 

Our next result gives the asymptotic variance of (/, //^) for / e B{W''). 

Theorem 2.2 Suppose ^ is k{'x.) -homogeneously stabilizing for K-almost all x e R'^. 
Suppose also that ^ satisfies the moments conditions (2.2) and (2.3) for some p > 2, 
and is power-law stabilizing for k of order q for some q with q > p/{p — 2). Let 
f G B{M.'^), and assume that Assumption Al or A3 holds, and that Assumption A4 
or A5 holds. Then 

lim A-Var[(/,4)] = / /(x) V«(x, «:(x))«:(x)dx, (2.14) 
with V^{x, a) given by V^{x., 0) = 0, and for a > by 

+a f (E[eS;^(Hr,M'')eS;^'(-z + ^a'^,K'^)]-(E[e^;^(^a,M'^)])2)rfz. (2.15) 

Also, the integral in (2.15) converges for K-almost all x, and the right hand side of 
(2.14) is finite. 

Our next result is a central limit theorem for the random field (A~^/^(/,7I^), / e 
B{W^)), or (A-V2^/,74),/ e -B(R'^)). We list some further assumptions. 

A6: For some p > 2, ^ satisfies the moments conditions (2.2) and (2.3) and is 
exponentially stabilizing for k. 

A7: For some p > 3, satisfies the moments conditions (2.2) and (2.3) and is 
power-law stabilizing for k of order q for some q > d{150 -\- 6/p). 

Theorem 2.3 Suppose \\k.\\oo < 00 and k, has bounded support. Suppose ^ is 
nip^) — homogeneously stabilizing at x for K-almost all x e , satisfies either AA 
or A5, and satisfies either A6 or A7. 

Then the finite- dimensional distributions of the random field 
(A~-^/^(/, /i^), / G B(R'^)) converge weakly as \ ^ 00 to those of a mean- zero finitely 
additive Gaussian field with covariances given by Jj^^ /i(x)/2(x)V^(x, fi;(x))fi;(x)(ix, 
with V^{x,a) given by (2.15). 

If also Al holds, the finite- dimensional distributions of the random field 
(A~^/^(/, 7l|), / G B(W^)) converge weakly as X 00 to those of a mean- zero finitely 
additive Gaussian field with covariances given by Jjjd /i(x)/2(x)l/^(x, «;(x))/«(x)(ix. 
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The corresponding results for the random measures lA^ ^ require some further condi- 
tions. These extend the previous stabihzation and moments conditions to binomial 
point processes. Our extra moments condition is 

inf sup sup E [^x(^,T;A:mUA*,M.'^Y] < oo, (2.16) 

X>l,xeR'',AeS3 (l-£)A<m<(H-£)A 

We give strengthened versions of A6 and A7 above, to include condition (2.16) 
and binomial stabilization. 

A6': For some p > 2, ^ satisfies the moments conditions (2.2), (2.3), (2.16), and is 
exponentially stabilizing for k and binomially exponentially stabilizing for k. 

AT: For somep > 3, ^ satisfies the moments conditions (2.2), (2.3) and (2.16), and 
is power-law stabilizing and binomially power-law stabilizing for k of order q 
for some q > d(150 -|- 6/p). 

For X e IR'^ and a > 0, set 

a) E [^^'CHa, R')] +a [ E [C^^{ni'^', R') - C^^iHa, R')]dy. (2.17) 

Theorem 2.4 Suppose \\k\\oo < oo and k has bounded support. Suppose ^ is 
k{'k) — homogeneously stabilizing at x for K-almost all x e M"^, satisfies Asuumption 
A4: or A5, and also satisfies A6' or AT. 

Then for any sequence (A(n),n e N) taking values in (0, oo), such that 
limsup„_^^ n~-'^/^|A(n) — n\ < oo, we have for f e S(R'^) that 

\unn-We.r{f,4 )= f /(x) V«(x, «(x))«(x)dx 

-(^^/(x)(5(x,/^(x))/^(x)dx^ , (2.18) 

and the finite- dimensional distributions of the random field (^^^^^(/i „)i / £ 

B{M'^)) converge weakly as n ^ oo to those of a mean- zero finitely additive Gaussian 
field with covariances given by 

^/i(x)/2(x)y^(x,«(x))«;(x)dx 
- ^ /i(x)5(x, «(x))«(x)dx^ /2(y)5(y, «(y))«(y)dy (2.19) 
with V^{'x.,X) given by (2.15). 
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If in addition, assumption Al holds, then (2.18) holds for f G B{W^) and the 
finite- distributions of the the random field {X{n)~^^'^ {f ,v\^^-^ _J, f € B{W^)) converge 
weakly as n ^ oo to those of a mean-zero finitely additive Gaussian field with 
covariances given by the expression in (2.19). 

Remarks. Theorems 2.2, 2.3 and 2.4 resemble the main results of Baryshnikov 

and Yukich [3], in that they provide central limit theorems for random measures 
under general stabilization conditions. We indicate here some of the ways in which 
our results extend those in [3]. 

In [3], attention is restricted to cases where assumption Al holds, i.e., where the 
contribution from each point to the random measures is a point mass at that point. 
It is often natural to drop this restriction, for example when considering the volume 
or surface measure associated with a germ-grain model, examples we shall consider 
in detail in Section 7.2. 

Another difference is that under Al, we consider bounded test functions in B{M.'^) 
whereas in [3], attention is restricted to continuous bounded test functions. By 
taking test functions which are indicator functions of arbitrary Borel sets Ai, . . . , Am 
in Mf^, we see from Theorem 2.3 that under Assumption Al, the joint distribution of 
(A~^/^/2^(Aj), 1 < i < m) converges to a multivariate normal with covariances given 
by jj^.^js^. V^^(fi;(x))/t(x)(ix, and likewise for i^x(n)n Theorem 2.4. This desirable 
conclusion is not achieved from the results of [3] , because indicator functions of Borel 
sets are not continuous. When our assumption Al fails, for the central limit theorems 
we restrict attention to almost everywhere continuous test functions, which means 
we can still obtain the above conclusion provided the sets A^ have Lebesgue-nuU 
boundary. 

The de-Poissonization argument in [3] requires finiteness of what might be called 
the radius of eodernal stabilization; see Definition 2.3 of [3]. Loosely speaking, an 
inserted point at x is not affected by and does not affect points at a distance beyond 
the radius of external stabilization; in contrast an inserted point at x is unaffected 
by points at a distance beyond the radius of stabilization, but might affect other 
points beyond that distance. Our approach does not require external stabilization, 
which brings some examples within the scope of our results that do not appear to 
be covered by the results of [3]. See the example of germ-grain models, considered 
in Section 7.2. 

In the non-translation-invariant case, we require ^ to have almost everywhere 
continuous total measure, whereas in [3] the functional ^ is required to be in a class 
SV(4/3) of 'slowly varying' functionals. The almost everywhere continuity condition 
on ^ is usually easier to check. 

We assume the underlying density function k is almost everywhere continuous, 
and for Theorems 2.3 and 2.4 but not for Theorem 2.2, we assume it has compact 
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support. In contrast, in [3] it is assumed that k has compact convex support and is 
continuous on its support (see the remarks just before Lemma 4.2 of [3]). 

Our moments condition (2.3) is simpler than the corresponding condition in [3] 
(eqn (2.2) of [3]). Using A7 and A7' in Theorems 2.3 and 2.4, we obtain Gaussian 
hmits for random fields under polynomial stabilization of sufficiently high order; the 
corresponding results in [3] need exponential stabilization. 

We spell out the statement and proof of Theorems 2.3 and 2.4 for the setting of 
marked point processes setting (i.e. point processes in M*^ x rather than in W^), 
whereas the proofs in earlier works [3, 16] are given for the setting of unmarked point 
process (i.e., point processes in R''). The marked point process setting includes many 
interesting examples such as germ-grain models and on-line packing, and generalizes 
the unmarked point process setting because we can always take Ai to have a single 
element and then identify R'^ x with R'^, to recover results for unmarked point 
processes from the general results for marked point processes. 

Other papers concerned with central limit theorems for random measures include 
Henrich and Molchanov [9] and Penrose [14] . The setup of [9] is somewhat different 
from ours; the emphasis there is on measures associated with germ-grain models 
and the method for defining the measures from the marked point sets (cqns (3.7) 
and (3.8) of [9]) is more prescriptive than that used here. In [9] the underlying point 
processes are taken to be stationary point processes satisfying a mixing condition 
and no notion of stabilization is used, whereas we restrict attention to Poisson or 
binomial point processes but do not require any spatial homogeneity. 

The setup in [14] is closer to that used here (although the proof of central limit 
theorems is difi:erent) but has the following notable differences. The point pro- 
cesses considered in [14] are assumed to have constant intensity on their support. 
The notion of stabilization used in [14] is a form of eoctemal stabilization. For the 
multivariate central limit theorems in [14] to be applicable, the radius of external 
stabilization needs to be almost surely finite but, unlike in the present work, no 
bounds on the tail of this radius of stabihzation are required. The test functions in 
[14] lie in a subclass of S(R'^), not i?(R'^). The description of the hmiting variances 
in [14] is different from that given here. 

In most examples, the sets F;^ are all the same as F. However, there are cases 
where moments conditions such as (2.7) and (2.16) hold for a sequence of sets F;^ 
but would not hold if we were to take F;s^ = F for all A. See, e.g. [15]. 

Last but not least, we note that our results carry over to the case where ^(x, t; X, •) 
is a signed measure with finite total variation. The conditions for the theorems re- 
main unchanged if we take signed measures, except that if ^ is a signed measure, the 
moments conditions (2.2), (2.3), (2.7) and (2.16) need to hold for both the positive 
and the negative part of ^. The proofs need only minor modifications to take signed 
measures into account. 
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3 Weak convergence: the objective method 



In the proof of Theorems 2.1 2.2, and 2.4, most of the work required is to prove 
convergence of first and second moments. A step in this direction is to obtain certain 
weak convergence results, namely Lemmas 3.3, 3.5, 3.6, 3.7 and 3.8 below. It is 
noteworthy that in all of these lemmas, the stabilization conditions used always refer 
to homogeneous Poisson processes on R'^; the notion of exponential stabilitization 
with respect to a non-homogeneous point process is not used until later on. 

To prove these lemmas, we shall use a version of what is sometimes called the 
'objective method' [1, 22], whereby convergence in distribution (denoted — >) for 
a functional defined on a sequence of finite probabilistic objects (in this case, re- 
scalcd marked point processes), is established by showing that these probabilistic 
objects themselves converge in distribution to an infinite probabilistic object (in this 
case, a homogeneous marked Poisson process), and that the functional of interest 
is continuous. We can then use the Continuous Mapping Theorem ([6], Chapter 1, 
Theorem 5.1), which says that if /i is a mapping from a metric space E to another 
metric space E', and Xn are E'-valued random variables converging in distribution 
to X which lies almost surely at a continuity point of h, then h{Xn) converges in 
distribution to h{X). 

A point process inW^ x M. is an £-valued random variable, where C denotes the 
space of locally finite subsets oiW^xM. We use the following metric on C: 

D{A, A') = {max{K e N : An [Bk x M] = A' n [Bk x M]})~^ . (3.1) 

With this metric, £ is a metric space which is complete but not separable. In the 
unmarked case where }A has a single element, our choice of metric is not the same 
as the metric used in Section 5.3 of [22]. Indeed, for one-point unmarked sets our 
metric generates the discrete topology rather than the Euchdean topology. 

To prove the weak convergence of point processes, we use a refinement of a 
coupling method used in [18]. In particular, we shall use a device which we here call 
the pivoted coupling. 

Let xq G and let A > 0. Let Ti"*" denote a homogeneous Poisson process 
of unit intensity in x M. x [0, oo). Let 7^^ denote the image of the restriction 
of 7i+ to the set {(x, t, s) E x M. x [0, cxd) : s < A«(x)}, under the mapping 
(x, t, s) f— > (x, t). For a > 0, let 7i^, denote the image of the restriction of H.^ to the 
set {(x, t,s) X M. X [0, oo) : s < Aa)}, under the mapping 

(x,i,s)^(A^/'^(x-xo),s). 

Then by the Mapping Theorem [11], V'x has the same distribution as V\ while Ti.'^ 
has the same distribution as Tia- We shall refer to V'x and as coupled realizations 
of V\ and Tia with pivot point xq. 
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Lemma 3.1 Suppose x e R'^ is a continuity point of k, and suppose (y(A), A > 0) 

is an -valued function which tends to x as A — oo. // and '^^^(x) ^'^^ 
coupled realizations ofV\ and Tin^x.) with pivot point y{\), then for any K e (0, oo), 
as A ^ oo, 

P[AV'^(-y(A) + V'^) n {Bk xM)^ K(x) n {Bk x M)] ^ 1. (3.2) 
Proof. The number of points of the point set 

{\"\-y{\) + AK(x)) n {Bk x M) 

equals the number of ponts (X, T, S) of Ti^ with X G -B;^-i/d^(y(A)) and with either 
Afi;(x) < 5* < Afi;(X) or Afi;(X) < 5* < Afi;(x). This is Poisson distributed with mean 

A / |k(z) — K(x)|dz 

which tends to zero by the assumed continuity of k{-) at x. □ 

In second moment computations, we are interested weak convergence, not only 
of a point process in £, but also of a pair of point processes in £ x £. The limiting 
object in this case will be a pair of independent homogeneous Poisson processes, 
and we need notation for this. To this end, for a > and 6 > let TY^ denote a 
homogeneous Poisson process on M'^ x TVf, independent of 7^^. 

Our first weak convergence lemma is concerned with the point processes X^- In 
this result, we assume these point processes are coupled togeteher in the natural 
way (this was not needed for the statement of results in Section 2). To do this, we 
let (Xi, Ti), (X2, T2), . . . denote a sequence of independent identically distributed 
random elements oiW^xM. with distribution k x , and assume the point processes 
'^m, Tn>l are given by 

— {(^1) Ti), (X2, T2), . . . , (X^, Tm)} (3.3) 

for each m. In proving the lemma, we use the notation = for equality of distribution. 

Lemma 3.2 Suppose (x, y) G M'^xM'^ with /t(x) > 0, /t(y) > 0, with «(■) continuous 
a^ X and at y, and with x 7^ y. Let (A(/c), i{k),m{k))k(z^ be a ((0, 00) x'NxN)-valued 
sequence satisfying X{k) — > 00, and £{k)/X{k) — > 1 and m{k)/\{k) — > 1 as A; —> 00. 
Then as k ^ 00, 

{\{kf/\-^ + A',(fc)), \{kfl\-^ + A'^(,)), \{kf'\-y + Xmi^k)), 
\{kfl\-^ + XJ^^;), \{kf/\-^ + A'^'J), \{kf'\-y + A'^'J))) 

— CHk(x), 'Hk(x), '^^(y), ?iK(x), 'Hk(x), '^^(y)). (3.4) 
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Proof. In this proof we ease notation by suppressing all mention of the marks in 
the notation, and write for example X for (X,T). Moreover, we suppress mention 
of the parameter k and write simply A for \{k), i for i{k), and m for m{k). 

We use the following coupling. Suppose we are given A. On a suitable probability 
space, let V and V be independent Poisson processes on R*^ with intensity function 
with each point carrying an independent Al -valued mark with distribution 
Hm- Let Yi,Y2, ... be independent random d-vectors with distribution k in 
carrying marks with distribution ^j^, independent of V and V. 

Let V' be the point process consisting of those (marked) points of V which lie 
closer to x than to y (in the Euclidean norm), together with those (marked) points 
of V which he closer to y than to x. Clearly V' is a (marked) Poisson process of 
intensity Ak(-) on W^. 

Let '^^'^(x) ^"^^ ^K(y) (marked) homogeneous Poisson point processes in M'^, of 
intensity /t(x) and /t(y) respectively. Assume '^^^{x) ^"^^ '^K(y) independent of 
each other and of (Yi, Y2, Y3, Y3, . . .). Assume also that '^^^(x) is coupled to V by 
the pivoted coupling with pivot at x (see Lemma 3.1), and 'H'^^^^ is coupled to V by 
the pivoted coupling with pivot at y. 

Let N denote the number of points of 7^' (a Poisson variable with mean A). 
Choose an ordering on the points of V' , uniformly at random from all N\ possible 
such orderings. Use this ordering to list the points of V as Wi, W2, . . . , Wjv- Also, 
set Wat+i = Yi, Wiv+2 = Y2, Wiv+3 = Y3 and so on. Set 

X', := {Wi, . . . , W,}, - {Wi, . . . , W^}. 

Clearly {X^X:^) ^ {X,,X„,), and (K(x),K(y)) = (H.(x), H.(y)) 

Let G N, and let 5 > 0. Let denote the volume of the unit ball in d 
dimensions. Define the events 

E := {Xl^ n 5A-v.i,(x) = n 5,-i/.^(x)}; 
F := + V)) HBk^ K(x) n Bk}. 

Event E occurs unless, either, one or more of the (N — m)'^ "discarded" points of V', 

or, one or more of the (m — iV)"*" "added" points of {Yi, Y2, . . .} lies in B)^-i/dj^{x^) . 
For each added or discarded point, the probability of lying in Bx-i/dj^{x) is at most 
6'||fi;||oo-^'^/A. Thus, for k large enough so that |m — A| < 5X, we have 

P[E^] < P[\N - A| > 5X] + {25X)e\\K\\^K^/X 

which is less than 3(56'||fi;||oo-f^'^ for large enough k. Hence, PlE'^] — >• as /c — > 00. 
Moroever, by (3.2) we also have P[F'^] — > as /c — > 00. 
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Assuming A is so large that |x — y| > 2 A ^I'^K, ii E f] F occurs then 

K(x) n 5x = (Ai/'^(-x + V)) n Bk 
= A^/'^((-x + v)n B^-./aj,) = +{vn B,-i/.^(x))) 

= AV'^(-x+(A';ni?,-v.i,(x))) 

= A^/'^((-x + A-;) n = (A^/'^(-x + A-J) n Bk 

so that -D('H^(x), AV''(-x + ^"4)) < Hence, for any K we have 

P[D(K(x),A^/'^(-x + AJ) > l/ii-] ^ 0, 

Similarly, we have 

max{P[Z}(K(x),A^/''(-x + a;)) > 1/K], P[DiK^^^, X'/'i-y + X;,,)) > 1/K], 

Combining these, we have the required convergence in distribution. □ 

For subsequent results, it is useful to define the region 

To := {x G int(r) : fi;(x) > 0, «;(■) continuous at x} (3.5) 

Lemma 3.3 Suppose (x, y) G Fq x Fq, with x 7^ y. Suppose also that R{x, T; TCk{x.)) 
and R{y,T;HK{y)) are almost surely finite. Suppose {\{m))m>i is a {0, 00) xN- valued 
sequence with \{m)/m — > 1 as m — >• 00. Then for Borel A C as m ^ 00 we 
have 

CamIx, T- X + X-^I^A) ^ ^<'{n,^^^,A), (3.6) 

and 

(6m(x, T; Xl'^', X + A-^Z-^A), eAM(y, T'; X^^ , y + X-''<'A)) 

^ (e^^^(^«(x), A), e^^'(7Y«(y), ^)). (3.7) 

Proof. Given A, define the mapping /i^ x : x £ — > R and the mapping h\ : 
(M X £ X M X £) ^ ]R2 by 

Vx(t^) =e(x,t;x + A,x + A); (3.8) 
h\{t,X,t',X') = {hA,^{t,X)MAt\^'))- 
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Since i?(x, T; 7Yk(x)) < oo a.s., the pair (T, 7Y«;(x)) lies a.s. at a continuity point 

of ^'A,x! where the topology on x >C is the product of the discrete topology on 
Ai and the topology induced by our metric D on C, defined at (3.1). Similarly, 
(T, 7iK{^),T', 7i«;(y)) lies a.s. at a continuity point of h\. We have by definition of 
that 

a(x, T; Xm, X + A-i/'^A) = HaAT, X'/'i-^ + X^)); 

= hl{T, aV'^i-x + xy'^'),T', xn-y + x^:^))- 

By Lemma 3.2, we have (T, A^/'^(-x + A"^)) ^ (T, H«(x)) so that (3.6) follows by 
the Continuous Mapping Theorem. Also, by Lemma 3.2, 

(T, X'/\-^ + A"^'^'), T', AV'^(-y + X^^')) ^ (T, H.(x), T', 7i«(y)) 
so that (3.7) also follows by the Continuous Mapping Theorem. □ 

The next lemma is one of several with the purpose of comparing (with reference 
to a test function /) the measure ^a(x; X^ ■) to the corresponding point measure 
^^(x; X ^ ■) so as to derive results when Assumption A2 or A3, rather than Al, is the 
case. In proving such results, we repeatedly use the notation 

0,(x) sup{|/(y) - /(x)| : y e 5,(x)}, for e > 0, (3.9) 

and for / e B{W^) we write ||/||oo for sup{|/(x)| : x e W^}. Also, recall (see e.g. 
[13], [20]) that X e M'^ is a Lebesgue point of / if s''^ J^^^^^ |/(y) - /(x)|dy tends 
to zero as £ J, 0, and that the Lebesgue Density Theorem tells us that almost every 
X e R'' is a Lebesgue point of /. 

Lemma 3.4 Let x e Fq, and suppose that R{'x,T;TIk{^)) < oo almost surely. Let 
y e R*^ with y 7^ X. Suppose that f G B(M'^), and suppose either that f is contin- 
uous at X, or that Assumption A2 holds and x a Lebesgue point of f . Suppose 
{X{m))m>i is a (0, oo) x N-valued sequence with X{m)/m — >• 1 as m ^ oo. Then as 
m — > OO; 

(/, eAM(x, T- X^) - ei(^)(x, T; A"^)) (3.10) 

and 

if, Cxim)i^, T- X^^') - CximM, T; Xl'^')) (3.11) 



17 



Proof. In this proof, we suppress the mark in the notation, writing x for (x, T), and 
,^(x; X, A) for ^(x, T; X, A), and so on. Also, we write A for A(m). The left side of 
(3.10) is equal to 



/ 



(/(z)-/(x))eA(x;A'^,rfz). 



(3.12) 



Given X > 0, we split the region of integration in (3.12) into the complementary 
regions S;^-i/ci;^(x) and R*^ \ S;^-i/d;^(x). Consider the latter region first. By (3.6) 
we have 



\B 



(x) 



(/(z)-/(x))a(x;A'^,dz) 



<2||/|UeA(x;A'^,M'^\B,-V.^(x)) 



where the limit is almost surely finite and converges in probability to zero as X — > oo. 
Hence for £ > 0, we have 



lim lim sup P 



(/(z)-/(x))eA(x;A'^,dz) 



> e 



= 0. (3.13) 



Turning to the integral over 5;^-i/d^(x), we consider separately the case where 
/ is continuous at x, and the case where A2 holds and x is a Lebesgue point of /. 
To deal with the first of these cases, observe that 



(/(z)-/(x))eA(x;A'^,c?z) 



<<Px-i/.K{^)U^;Xm,R'') (3.14) 



and if / is continuous at x, then 0;^-i/d^(x) 0, while ^x{-x, Xm,M.'^) converges in 
distribution to the finite random variable ^^^{7iK{x),^'^) by (3.6), and hence the 
right hand side of (3.14) tends to zero in probability as m — > oo. Combined with 
(3.13), this gives us (3.10) in the case where / is continuous at x. 

Under Assumption A2, for Borel A C M*^, the change of variables z = x + 
A-Vd(y _ x) yields 



eA(x; X,A)= f e'(x; X + X'/'{-^ + X), y)dy 

^x+Ai/d(-x+A) 

^xj C'(x;x + A^/'^(-x + A'),x + A^/'^(z-x))dz. 
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Hence, under A2, 



= A 



(/(z)-/(x))eA(x;AW,o?z) 

/ (/(z) - /(x))e'(x; X + AV<^(-x + A'^), X + AV'^(z - x))dz 

<XoA / |/(z)-/(x)|dz 



and if additionally x is a Lebesgue point of / then this tends to zero. Combined 
with (3.13), this gives us (3.10) in the case where A2 holds and x is a Lebesgue 
point of /. 

The proof of (3.11) is similar; we use (3.7) instead of (3.6). □ 

By combining Lemmas 3.3 and 3.4, we obtain the following, which is the main 
ingredient in our proof of the Law of Large Numbers in Theorem 2.1. 

Lemma 3.5 Suppose (x, y) G FoxFo, with x 7^ y. Suppose also that i?(x, T; T-Ik{x)) < 
00 and i?(y, T; 7YK(y)) < 00, almost surely. Let f G B(M.'^) and suppose that either 
Al holds, or A2 holds and k is a Lebesgue point of f, or 'x. is a continuity point of f. 
Suppose {X{m))m>i is a (0, 00) x N-valued sequence with \{m)/m — > 1 as m ^ 00. 
Then as m ^ 00, 



if, eAM(x, T; Xrn)) ^ /(x)ej;^(7i,(^), R'^) 



(3.15) 



and 



(/, 6(„o (x, T; X^'^')){f, 6(H (y, T'; X:^' )) 



V 



(3.16) 



Proof. Note first that 

(/, a(^)(x, t- X)) = /(x)^,(^)(x, t- X, W). 

Hence, in the case where Al holds (i.e., ^ = ^*), (3.15) is immediate from the case 
= IR'^ of (3.6), and similarly (3.16) is immediate from (3.7). 
In the other two cases described, we have (3.10) by Lemma 3.4. Combining 
this with (3.15) for the case with ^ = we see by Slutsky's theorem (see, e.g., 
[13]) that (3.15) still holds in the other two cases. Similarly, since (3.16) holds when 
^ = by (3.11) and Slutsky's theorem we can obtain (3.16) in the other cases too. □ 

The next two lemmas are key ingredients in proving Theorem 2.2 on convergence 
of second moments. 
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Lemma 3.6 Suppose x e Fq. Suppose also that ^ is k{x) -homogeneously stabilizing 
at X, and that Assumption A4 or A5 holds. Then for any z e R*^, we have 

eA(x + A-i/'^z,T;7^A,M') ^e^;^(^«(x),M") as A ^ oo. (3.17) 
Also, if f & B^U!^) and f is continuous at x, then 

{f,U^ + \-^/\T-Vx)) ^ /(x)e^^(H«(x),M'^) as A ^ oo. (3.18) 

Proof. The proof of (3.17) is related to that of Lemma 3.3. Given x and z, set 
Va := X + A-V'^z. For Borel A C R'^, define : R'^ x Al x £ ^ R by 

g(A(w, t, X) = ^(x + w, t; X + w + A", X + w + A). 

Then 

a(vA, T; Vx, VA + A-^Z-^A) = gA{X-'^\ T, \"\-vx + Pa))- 

Taking our topology on R'^ x x £ to be the product of the Euclidean topology 
on R^, the discrete topology on M. and the topology induced by the metric D on L 
which was defined at (3.1), we assert that as A ^ cxd, 

(A-V<iz,r,AV<i(-VA + n)) ^ (0,r,7i,(x)). (3.19) 

To see this, for each A, let and '^^^(x) coupled realisations of V\ and 

obtained by the pivoted coupling with pivot at va. Then for £ > 0, by Lemma 3.1 
we have P[D{V')^,n'^^^^) > s] ^ as A ^ oo. This gives us (3.19). 

If Assumption A4 (translation invariance) holds, then the functional fi'yi(w, t, X) 
does not depend on w, so that 5'a(w, t, X) = 5'a(0, t, X) and by the assumption that 
^ is «;(x) -homogeneously stabilizing at x, we have that (0, T, 'Hk(x)) almost surely 
lies at a continuity point of the functional gA- 

If, instead. Assumption A5 (continuity) holds, take ^4 = R*^ or ^4 = Bk or 
A = W^\Bk, with K > Ki and given in Definition 2.5. Then by the assumption 
that C, is /t(x) -homogeneously stabilizing at x (see (2.10)), with probability 1 there 
exists a finite (random) r] > such that for D{X ,7iK{:x.)) < ^, and for |w| < r], 

5a (w, T, A") = e(x + w, T; X + w + (7i„(,) n Si/^), x + w + A) 
^ e(x, T; X + (n^^^) n Si/^), X + A) = 5a(0, T, 7i«(,)) as w ^ 0. 

Hence, (0, T, 7^fj(x)) almost surely lies at a continuity point of the mapping gA in 
this case too. 
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Thus, if A is R"^ or Bk or IR*^ \ Bk, for any K under A4 and for K > Ki under 
A5, the mapping qa satisfies the conditions for the Continuous Mapping Theorem, 
and this with (3.19) gives us 



Uvx, T; Vx, va + X-'/'A) ^^^(n,^^),A) as A ^ oo 



(3.20) 



Taking A = M'^ in (3.20) gives us (3.17). 

Now suppose that / is continuous at x. To derive (3.18) in this case, note first 
that 

{f,U^x,T;Vx)-Cxi^x,T;V,))= [ (/(w) - /(v,))e,(v,, T; P,, dw). 
Given K > 0, by (3.20) we have 



{f{w)-f{^^x))U^x,T■,Vx,dw) 



;,-l/dx(VA) 

< 2||/||ooeA(vA, T; Vx, M'^ \ S;,-v.^(va)) 



V 



where the hmit is almost surely finite and converges in probability to zero as X — > oo. 
Hence for £ > 0. we have 



lim lim sup P 

K-*oo x-^oo 



/Md\B^_l/d^(vA) 

Also, given > 0, it is the case that 



{f{w)-f{vx))U^x,T;Vx,dw) 



> e 



0. (3.21) 



B;,-l/d^(VA) 
< 



{f{w)-f{vx))U^x,T;Vx,dw) 
20A-v<^(K+|z|)(x)6(vA,r;PA,R') 



(3.22) 



and by continuity of /, 0A-i/d(x+|z|)(x) ~^ while ^xi'^x^T^Vx^R^) converges in 
distribution to the finite random variable $,^^{Ti.K(^y_),W^) by (3.17), and hence the 
right hand side of (3.22) tends to zero in probability as A ^ oo. Combined with 
(3.21), this gives us 



(/, 6(VA, T; Vx) - Cxi^x, T- Vx)) 0. 
Also, by (3.17) and continuity of / at x, we have 

(/,a(vA,r;PA)) ^ /(x)e^^(H«(x),M'^), 
and combined with (3.23) this yields (3.18). 



(3.23) 



□ 
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Lemma 3.7 Suppose ^ satisfies Assumption AA or A5. Then for Lehesgue- almost 
all X e Fq and all z e IR*^, as \ ^ oo we have 



^ eS;^(K£>')e^;^'(-z + nZ,X). (3.24) 



Also, for f e B{M.'^), if f is continuous at x, then 

if, 6(x, T; VT'-'"^'^')) X (/, + A-V'^z, T'; P.^'^)) 



^ + nZ,X). (3.25) 



Proof. Again write y\ for x + A ^/"^z. Let ^4 C M*^ be a Borel set. Define the function 
: M'' X A< X X £ ^ by 

^a(w, t, t', X) = (e(x, x + A'^'*', X + ^), 
^(x + w, i'; X + w - z + X^''^, X + w + A)). 

Then 

(6(x, T; PJ-^', X + A-V-^A), eA(vA, T'- V^'^, v, + A'^-^^)) 
= (e(x, T; X + A^/'^(-x + PJ-^'), X + A), 

e(v,, T'; VA + A^/'^(-x - \-'/'z + Vt'^), + A)) 
= 5A(A-V'^z,r,r',AV'^(-x + P,)). 

Under A5, let us restrict attention to the case where A is W^, Bk or W'' \ Bk with 
K > K\. Then under either A4 or A5, by similar arguments to those used in proving 
Lemma 3.6, (0, T, T', 'Hk{^)) lies almost surely at a continuity point of qa, and since 

A~^/'^(— X + Pa) "^/{(x), the Continuous Mapping Theorem gives us 



(a(x, T; PJ-^', X + A-i/'^A), eA(vA, T'; P^^, va + X-'/'A)) 



^^(0, T, T', H«(x)) = (^^^(K£U)> + ^)) (3-26) 



as A ^ oo. Taking ^ = R'^ gives us (3.24). 

Now suppose / is continuous at x. To prove (3.25) in this case, observe that for 
X > 0, we have 

|(/,6(x,T;Pj-^') -a(x,T;Pj-^'))| < 0A-v.;,(x)a(x, T; Pj-^', M'^) 

+2||/|UeA(x,r;Pj-^',M'^\S,-v.x(x)) (3.27) 
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The first term in the right hand side of (3.27) tends to zero in probabihty for any 
fixed K, by (3.26) and the fact that is almost surely finite. Also by 

(3.26), the second term in the right hand side of (3.27) converges in distribution, 

as A ^ oo, to 2II/II c^^T(^^,T' 

K ^ oo. Hence, by (3.27) we obtain 



^^^iK^y^'' \ Bk), which tends to zero in probability as 



(/,eA(x,T;Pj-^')-a(x,T;P,---)) 



VA,T'x 



0. 



(3.28) 



We also have 



< 



(/(y)-/(x))6(vA,r';P,-'^ciy) 



+2\\f\U,{^,,T';V: 



\5^-i/d^(vA)). 



(3.29) 



By (3.26) and the assumed continuity of / at x, the first term in the right side of 
(3.29) tends to zero in probability for any fixed K, while the second term converges 
in distribution to 2 1 1 / 1 1 ' - ■ -^O'^ 
as X — > oo. Hence, as A ^ oo we have 



zero in probability 



(/, 6(VA, T'- Vr ) - a(vA, T'- Vr )) 0. 
By continuity of / at x, and the case A = R'^ of (3.26), we have 

((/, a(x, T- Vl^'^')), (/, a(vA, T'; V-/))) 
Combining this with (3.28) and (3.30) yields (3.25). 



(3.30) 



□ 



The following lemma is a refinement of Lemma 3.5 and is proved in the same 
manner as that result. It will be used for dc-Poissonizing our central limit theorems. 
To ease notation, we do not mention the marks in the notation for the statement 
and proof of this result. 

Lemma 3.8 Let (x, y) G Fq with x 7^ y, and let (z, w) G (R'^)^. Suppose either that 
Assumption Al holds, or that x and y are continuity points of f . Given integer- 
valued functions (£(A),A > 1) and (m(A),A > 1) with i{X) ~ A and m(A) ^ X as 
A — > OO; we have convergence in joint distribution, as A — > 00, of the 11 -dimensional 
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random vector 

((/,eA(x; A',)), (/,a(x; A'/)), (/, ^(x; A'^""''^)), {f,U^;Xr'""'^ U {y})), 

{f,U^;Xj), (/,eA(x; A"^)), {f,U^;X^, U {x + A-V'^z})), (/,eA(y; ^-J), 

(/, 6(y; A"^)), (/, ^,(y; A"^ U {x + X-'/'z})), 
if, 6(y; u {x + x-'/\ y + A-V'^w}))) 

to 

M-^), /(x)e^(H,(.), M-^), /(x)C(K(x), K'), /(x)C(K(x), K'), 

/(y)^^(^«(y), K^), /(y)e^(^«(y), K'), /(y)e^(^r(y), K')) • 

Proof. First, we assert that 

(eA(x; X,, R'),U^; A/, M'^), a(x; X^^''"''', R''),U^; Af +^"''^' U {y}, M'^), 

eA(x; A^, M'^), eA(x; A^, M"^), ^(x; A^ U {x + A-^-^z}, R'^), ^Aly; A'„„ R'^), 

eA(y; A'„^., R'), a(y; u {x + x-'/'z}, m'), 

Ut, U {x + A- V'^z, y + X-'/'wjX)) 

converges in distribution to 

(nniy) , R') , (7i.(y) , R'^) , (H^iy) , R'^) ) • 

This is deduced from Lemma 3.2 by a similar argument to the proof of Lemma 3.3. 
For example, considering just the third component, defining the mapping /i-Rd^x on 
M X Chj (3.8), we have 

eA(x; R<^) = h^.,^(T, AV'^(-x + 

= Vx(r,MuAV<i(-x + A',)), 

and by Lemma 3.2, (T, {z} U Ai/'^(-x + A^)) converges in distribution to (T, ) 
which is almost surely at a continuity point of h^d^^. Similar arguments apply for 
the other components and give us the assertion above. This assertion implies that 
the result holds under Al, i.e. when ^ — 
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Now let us drop Assumption Al, but assume that x and y are continuity points 
of /. Then by Lemma 3.4, 



(/, U^; X,) - e:(x; X,)) 0, (/, Xj) - ^^(x; Xj)) 0, 

(/, a(x; X^) - ei(x; X^)) 0, (/, 6(x; XI) - Cxi^; X.^)) ^ 0, 

(/, 6(y; - a(y; ^n.)) ^ o, (/, 6(y; x^^) - Cxiy; ^^)) o, 



and a similar argument to the proof of (3.28) (working with Xm instead of V\) yields 



(/, eA(x; Xr'-'^"^ U {y}) - ^^(x; xr'"'"^ U {y})) ^ 0, 
(/, eA(x; X--^'-'^'^ U {y}) - a(x; A'^'"''^ U {y})) ^ 0, 
(/, e,(y; A'^ U {x + A'^'^z}) - G(y; X-U{^ + A-^-^z})) ^ 0, 
(/, Ur, U {x + A-^/'^z, y + A-i/M) 
-a(y; u {x + A-^z-^z, y + A-i/'^w})) ^ 0. 



Combining these eleven convergence in probability statements with the fact that we 
have established our conclusion in the case where Assumption Al (^ = ^*) holds, 
and using Slutsky's theorem, we obtain our conclusion in the other case as well. □ 



4 Proof of Theorem 2.1 

Before proving Theorem 2.1, we give general expressions for the first two moments 
of (/, /i^) which we shall use again later on. By Palm theory for the Poisson process 
(e.g. a slight generahzation of Theorem 1.6 of [13]), we have 




)) 



p 



0. 



(4.1) 



and 




+2E J2 {f,U^,t;T^x)){f,Uy,u;'Px)) 



{(x,t),(y,«)}C7'A 



= \E[{fU^,T;Vx))'] 
+X'E [(/, U^, T; Vf'^')) if, U^', T'; ^))]. 



(4.2) 
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Combining (4.1) and (4.2), we have 



A-2Var(/, /ii) = \~'E[{f,U^,T-Vx))'] 
+E [(/, eA(X, T; Vf'^')){f, U^', T'; Vf'^))] 

-{E[{f,UX,T;rx))]y. (4.3) 

Also, by similar arguments, 

n-'E if, 4J = E (/, eA(X, T; (4.4) 

and 

n-^Var(/, 4 J = n'^E [(/, ^(X, T; A'„_i))2] 
+ (^^) E[(/,6(X,r; A-^^^f ))(/,eA(X',r'; A-^^J^))] 

-(E[(/,eA(x,r;A;_0)])'. (4.5) 

Recall that by definition, ^x{x,t; X,R'^) = for x G M'^ \ Tx, and (F^, A > 1) is a 
given nondecreasing family of Borel subsets of M.'^ with limit set F having Lebesgue- 
nuU boundary; in the simplest case Fa = M'' for all A. 

Proof of Theorem 2.1. First we prove (i) for the case q = 2. Assume that (2.2) and 
(2.3) hold for some p > 2. Let 'Hk{'k.) denote a Cox point process in x TVf, whose 
distribution, given X = x, is that of Hn[^). Set J := /(X)^^'^(7i^(x), M'')lr(X), 
and let J' be an independent copy of J. 

For any bounded continuous test function h on R, by (3.15) from Lemma 3.5, 
as A ^ cx) we have E [/i((/, ^a(X, T; Pa)))|X] E[J|X], almost surely Hence, 
E[/i((/,a(X,r;PA)))] ^E[/i(J)], so that 

(/,eA(X,r;PA)) ^ J. (4.6) 
Similarly, using (3.16) we obtain 

(/,eA(X,r;Pf '^'))(/,eA(X',T';Pf ^)) ^ J' J (4.7) 

Also, by (2.2) and (2.3) the variables in the left side of (4.6) and in the left hand 
side of (4.7) are uniformly integrable so we have convergence of means in both cases. 
Also, (2.2) shows that the first term in the right side of (4.3) tends to zero. Hence we 
find that the expression (4.3) tends to zero. Moreover, by (4.1) and the convergence 
of expectations corresponding to (4.6), A~^E (/, /i^) tends to E [J], and this gives us 
(2.5) for q — 2, under the assumptions of part (i) of Theorem 2.1. 
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Now consider the case q — 1. Assume (2.2) and (2.3) hold for some p > 1. We 

use a truncation argument. Define positive and negative parts and /~ of the 
test function / by /^(x) := max(/(x),0) and /^(x) := max(— /(x), 0). For K > 0, 
let truncated version of the measure defined by 

C^(x, t; X, A) := min(^(x, t] X, W^), K)1a{^). 

Then has total measure bounded by K. By the case q — 2 established above, 

A-^(r,/^r) ^ E [/+(X)(e^)^'^(7i.(x),M'^)lr(X)] . (4.8) 
Then (/+,r(x,i; AT)) = lim;^^^(/+, ^^(x, A")). Using (4.1), we have 

0<E[A-i(/+,/;)-A-^(/+,/if)] 
= E [(/+, a(X, T; Vx) - (X, T; P,))] 
< ||/||ooE[eA(X,T;P,,M'^)l{eA(X,T;P,,M'^) > i^}] 

which tends to zero a,sK ^ oo, uniformly in A, because the moments condition (2.2), 
p > 1, implies that the random variables ^xil^.T^Vx^W^) are uniformly integrable. 
Also, by monotone convergence, as K ^ oo the right side of (4.8) converges to 
E[/+(X)e^'^(7i,(x),M'')lr(X)]. Hence, taking X ^ oo in (4.8) yields 

A-^(r,/^f ) ^ E[/+(X)e'^(7i.(x),M'^)lr(X)], 

and a similar argument yields the equivalent statement with /+ replaced by /~. 
Combining these, and using linearity, we obtain 

A-i(/,/;) = A-i((/+,^f)-(r,^f)) 

^ E[(/+(X)e^'^(H«(x),M'^) - /-(X)e^'^(H.(x),M''))lr(X)] 

= E[e^'^(7i,(x),M'^)/(X)lr(X)]. (4.9) 

This gives us (2.5) for = 1 when Assumption Al hols. 

Now suppose A2 or A3 holds. By Palm theory, analogously to (4.1) we have 

EA-^ J]|(/,eA(x;PA)-a(x;PA))| 

= ^ «:(dx)E I (/, U^- Vx) - Vx))\. (4.10) 

and for almost every x e Fi, by (3.10) and (2.2), the integrand tends to zero and is 
bounded so that we have (2.6). Combining this with (4.9) gives us (2.5) for q — 1 
when Assumption A2 or A3 holds, completing the proof of part (i). 
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Next we turn to part (ii) with q — 2. Assume (2.7) holds for some p > 2. By 
Lemma 3.5, as n — > oo we have 

(/,eA(n)(X,T;A'„_0)^ J; (4.11) 
(/,6(n)(X,T; A'„^_:f ))(/,eA(„)(X',T'; A^^^J^)) ^ fj, (4.12) 

Also, by (2.7) and the Cauchy-Schwarz inequality, the variables in the left side of 
(4.11) and in the left side of (4.12) are uniformly integrable so we have convergence 
of means in both cases. Likewise, (2.7) shows that the first term in the right side of 
(4.5) tends to zero. Hence we find that the expression (4.5) (with A = A(n)) tends 
to zero. Also, by (4.4) and (4.11), n~'^'E {f,i/x(n),n) IE[J], and this gives us (2.8) 
for q = 2. 

The case g = 1 of part (ii) is deduced from the case g = 2 in the same manner 
as in part (i). □ 

5 Proof of Theorems 2.2 and 2.3 

Note that by definition, ^^(x, t; X, W^) = for x e R'' \ L. In the sequel, we fix a 
test function / e S(R'^). Set 



ax :=^E[(/,a(x,T;PA))']«^(x)cix (5.1) 

and 

/3a := / / (E [(/, eA(x, T; V^,^'-''"^''^')) (/, ^(x + X'/\ T'; P.^'^))] 
Jr JRd, 

-E [(/, eA(x, T; V,))]E [(/, eA(x + X-'/'z, T; Vx))]) 

xk(x)k(x + A~^/'^z)dzdx; (5.2) 

Lemma 5.1 For X> 0, it is the case that 

Var((/,4)) = A(aA + /9A). (5.3) 
Proof. The first term in the right hand side of (4.3) equals A^^cka. Thus, (4.3) yields 

A-Var((/,/.i)) -CA = X JJjE[{f,U^,T-,Vl'^')){f,UyX;V-x^))] 

-E [(/, eA(x, T; Vx))]E [(/, Uy, T; Vx))M^)K{y)dyd^, (5.4) 

and the change of variables y = x + A^^/'^z shows that this equals Px as given by 
(5.2). □ 
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Lemmas 3.6 and 3.7 establish hmits in distribution for the variables inside the 
expectations in the integrands in the expressions (5.1) and (5.2) for a\ and Px. To 
prove Theorem 2.2, we need to take these limits outside the expectations and also 
outside the integrals, which we shall do by a domination argument. It is in this step 
that we use the condition of stabilization with respect to non-homogeneous Poisson 
processes (Definition 2.4), via the following lemma, which is an estimate showing 
that the integrand in the definition (5.2) of /Sx is small for large |z|, uniformly in 
A. To ease notation, for x G M'^, z G and A > 0, we define random variables 
X = X^,^,x, Z = Z^,^,x, X' = X4 , , and Z' = by 

X := (/, 6(x, T; T^r'"'"''^')), Z := (/, 6(x + A'^Z-^z, T'; V^'"^)), (5.5) 
X' := {f,U^,T;Vx)), Z' - (/, eA(x + A-V<^z, T'; 7^,)). (5.6) 

Similarly, we define random variables X* — X'^^y^, Z* — Z^^^, X*' — X^^^ and 
^*' = <z,A,by 

X* :=eA(x,T;P^+^"'^'"'^',M'^), Z* := eA(x + A'^/^^z, T'; pf^, R'^), (5.7) 
X*' ■.= U^,T;Vx,m''), Z*' :=eA(x + A-i/'^z,T';PA,M'^). (5.8) 

We write a Ab for min(a, b) in the sequel. 

Lemma 5.2 Suppose that ^ satisfies (2.2) and (2.3) for some p > 2, and is power- 
law stabilizing for k of order q for some q > dp/{p — 2). Then there is a constant 
Ci, independent of \, such that for a// A > 1, x G supp(/€) and z G M.'^, 

|E [Xx,.,a^x,z,a] - E [X;,,JE < Ci(|z|-'^-(Vc^^) A 1); (5.9) 

|E K,.,a^x,z,a] - E K;,,,]E < Ci(|z|-<^-(Vt^i) A 1). (5.10) 

Proof. Let X := Xx,z,a and Z := Zx,z,a- Let X = ^l{ii^(x,r)<|z|/3} and let 
Z = ^l{K^(x+A-i/dz,T')<|z|/3}- Then X and Z are independent so that E [XZ] = 
E[X]E[Z], and 

E [XZ] = E [X]E [Z] + E [X{Z -Z)]+E [{X - X)Z] (5.11) 

while 

E [X']E [Z'] = E [X]E [Z] + E [X]E [Z' -Z]+E[X' - X]E [Z']. (5.12) 

By (2.3) and Holder's inequality, and the assumed power-law stabilization of order 

q > dp/{p — 2), there is a constant C2 such that 

E[|(X-X)Z|]=E[|XZ|l|«,(x)>|z|/3}] 
< {E[\X\P]y/P{E[\Z\P]y/P{P[Rx{^,T) > |z|/3])i-('/P) 

< Cadzl"'^"^'/^'^ A 1) (5.13) 
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and likewise 

E[\X{Z-Z)\] <C2(|z|-'^-(^/^2) Al). (5.14) 
By a similar argument using (2.2), there is a constant C3 such that 

max(E [\X' - X\],E [\Z' - Z\]) < C'sdzj-'^-^^/^^) ^ (5^5) 

Subtracting (5.12) from (5.11) and using (5.13), (5.14) and (5.15) along with the 
Cauchy-Schwarz inequality, we may deduce that there is a constant Ci, independent 
of A, such that for all A > 1, (5.9) holds. The argument for (5.10) is similar □ 

Proof of Theorem 2.2. Let / G B{W^). Let x G Fq (as defined at (3.5)). Assume 
also, either that Al holds, or that A3 holds and x is a continuity point of /; also 
assume A4 or A5 holds. By the case z = of (3.17) when Al holds, and the case 
z = of (3.18) when A3 holds, we have 

{f,U^,T;Vx)) ^ /(x)eS;^(7i,(.),M'^). (5.16) 

By (2.2), {(/, ^a(x; Pa))^ : A > 1} are uniformly integrable, and hence the conver- 
gence in distribution (5.16) extends to convergence of second moments to a limit 
which is bounded by (2.2) and Fatou's Lemma. Hence, by the dominated conver- 
gence theorem, ax given by (5.1) satisfies 

\imax= /(/(x)E[C^^(H,(x),M")])'/s;(x)rfx<oo. (5.17) 

Next we show convergence of the expression Px given by (5.2). To this end, set 

(7a(x, z) := (E [X.,.,a^x,z,a] - E [X^, K,,,x])<^ + A^'/'z); 
^^(x, z) := (E [X*^,^,xZ*.,.,x] - E [<.,a]E [<.,a])'^(x + A-^Z-^z). 

Suppose A3 holds. Then for almost all x e Fq and all z G M.^, with «;(x) > and k, 
continuous at x, by Lemma 3.7 we have as A ^ 00 that 

X.,.,aZ.,z,a f{^rCj{n%,R')^:f'{-z + K§yR'), (5.18) 

and the variables in the left side of (5.18) are uniformly integrable by (2.3), so that 
(5.18) extends to convergence of expectations. Likewise by (2.2) and (3.18), both 
E[X'^^^J and E[Z'^^^J converge to f {^)E [^^^ {ni^^yR% so we have under A3 
that 

lim (^a(x, z)) = ^oo(x, z), a.e. (x, z) G Fq x W^, (5.19) 

A— »oo 
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with 



-E[C^(H«(.),K'^)]^). 

By our assumptions on k{-) and F, 

/ «;(x)dx = 0. (5.20) 
ir\ro 

By Lemma 5.2, and the assumption that k is bounded, there is a constant C such 
that |5'a(x,z)| < Cdzj-'^-^i/^) A 1), for almost every (x, z) G Tq x M"^, and A > 1. 
Hence, by (5.2), (5.19), (5.20), and the dominated convergence theorem, we have 



Px— / 5'a(x, z)fi;(x)dzdx — > / / ^(^^(x, z)fi;(x)(izdx < oo 





Combined with (5.17) and (5.3), this gives us (2.14) under A3. 

Now suppose instead of A3 that Al holds. Then by (3.24), for almost all (x, z) e 
To X M*^ we have as A ^ oo that 

and the variables in the left side of (5.21) are uniformly integrable by (2.3), so that 
(5.21) extends to convergence of expectations. Likewise by (2.2) and (3.17), both 
E [X*'^ and E [Z^^'^,;,] converge to E [i'^{'Hl^^y M.% so we have under Al that 



lim(^*(x,z))=^:,(x,z), a.e. (x,z)eroX]R'^ (5.22) 



with 



^:,(x, z) «:(x)(E {-^ + n^^y^')] - E r^^in.^^'^X)?). 

By Lemma 5.2, the assumption that k is bounded, and (5.22) there is a constant C 
such that for almost every (x, z) with /t(x) > 0, 

|^*(x,z)| <C(|z|-<^-(Vf^) Al), l<A<oo. (5.23) 

If X e R'' is a Lebesgue point of /, then for any X > 0, by (5.10) we have 

hm / y:(x,z)(/(x + A-V''z)-/(x))dz = 0, 
Jbk 
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and combining this with (5.22) and the dominated convergence theorem gives us 



hm / ^*(x,z)/(x + A-V<^z)dz= / <7^(x,z)/(x)dz. (5.24) 
On the other hand, by (5.23) and the assumption that / is bounded, we have 
hm hmsup / |^*(x, z)/(x + - g*^{^, z)/(x)|dz = 

K-*oo x-^^ jRd\BK 

and combining this with (5.24) we have 

hm / <7l(x,z)/(x + A-V<iz)rfz= / gl{^,z)f{^)dz. 

By the Lebesgue density theorem, almost all every x e Fq is a Lebesgue point of /. 
Hence, under Al, by (5.2), (5.20), the dominated convergence theorem, 

/3a = / /(x) / /(x + A-^/'^z)y*(x,z)Ac(x)dz(ix^ / /(x)2 / g*^{^,z)K{^)dzd^, 

and combined with (5.17) and (5.3), this gives us (2.14) as required. □ 



For the proof of Theorem 2.3 (central limit theorem for random measures), we 
shall use results on normal approximation for (/, /I^), suitably scaled. In the case of 
point measures, these were proved by Stein's method in [19], and the method carries 
through to more general measures. Let $ denote the standard normal distribution 
function, and let AA(0, cr^) denote the normal distribution with mean and variance 
(T^ (if (j^ > 0) or the unit point mass at if cr^ = 0. 

Lemma 5.3 Suppose that n has bounded support and \\k\\oo < oo. Suppose that ^ is 
exponentially stabilizing and satisfies the moments condition (2.2) for some p > 2. 
Let f e B{W^), and q e (2, 3] with q < p- There exists a finite constant C depending 
on d, ^, K, q and f , such that for all X> 1 



sup 



P 



if, A) 



(Var(/,/.i))V2 



< t 



<C(logA)^<^A(Var(/,/ii)) 



(5.25) 



Proof. In the case where ^ = i.e. ^(x, A", •) is always a point mass at x, this 
result is Theorem 2.1 of [19]. If we do not make this assumption on ^, the proof in 
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[19] carries through with httle change, except that the Tx and T( of Section 4.3 of 
[19] should now be defined (following notation of [19]) by 

V(A) Ni 

^^ = EE(/'^a(a^'^«'^a)), 

V{X) 
i=l j=l 

Lemma 5.4 Suppose \\k\\oo < oo. Suppose for some p > 3 that ^ is power-law 

stabilizing of order q for some q > d{150 + 6/p) , and satisfies the moments condition 
(2.2). Let f G B(M.'^). Suppose that X^^VaT{f, fi\) converges, as X ^ oo, to a finite 
limit cr^. Then (/, A"^/^/!^) converges in distribution, as X ^ oo, to the J\f {0, a'^) 
distribution. 

Proof. In the case where ^ = this result is Theorem 2.2 of [19]. If we do not 
make this assumption on ^, the proof in [19] carries through with the same minor 
changes as indicated for Lemma 5.3 above. □ 

Proof of Theorem 2.3. Suppose ||/t||oo < oo and k has bounded support. Suppose 
^ is almost everywhere continuous, and is fi;(x)— homogeneously stabilizing at x for 
K-almost all x e IR*^. Suppose ^ satisfies either A4 or A5, and satisfies A6. 

Let / e S(R'^), and assume either Al or A3 holds. By Theorem 2.2, A-^Var(/, /x^) 
converges to the finite nonnegative limit Jp /(x)^V"^(x, K(x))/«((ix). If this limit is 
strictly positive, then the right hand side of (5.25) tends to zero, so that by Lemma 
5.3, A~^/^(/, 71^) is asymptotically centred normal with variance 
Jp /(x)^V^(x, K(x))K((ix). On the other hand, if this limiting variance is zero, then it 
is immediate from Chebyshev's inequality that A~^/^(/, converges in probability 
to zero. Hence for all / e B{R'^), or under Al for all / e B{R'^), we obtain 

^-'"{fA) ^AA(^0, j|^/(x)V«(x,Ac(x))/.(x)dx^ . (5.26) 

If A7 holds instead of A6, we obtain the same conclusion by using Theorem 2.2 and 
Lemma 5.4; note that in A7, since p > 3 the condition q > d{150 + 6/p) ensures 
that q > dp/{p — 2), so that Theorem 2.2 still applies here. 

Now consider an arbitrary finite collection of test functions fi, . . . , fj, each of 
them in i?(K'^). For arbitrary real constants bi,...,bj, application of (5.26) to 
/ = E/=i bjf) yields 

^6,.A-V2(/.,74) ^ AT 0,5^5^ / /,(x)/,(x)y^(x,«:(x))«(x)dx 
j=\ \ j=i k=i -^^ 
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so by the Cramer- Wold device (see e.g. [13], [6]), the variables \~^/'^{fj,Ji^^) , 1 < j < 
J, are asymptotically centred multivariate normal with covariance matrix having 
(j, A;)th entry Jp /j(x)/fe(x)y^(x, /t(x))/s;(x)(ix. This gives us the required conver- 
gence of random fields for fj G B{M.'^). If Al holds, the same argument gives the 
required convergence of random fields for fj G BiW^). □ 



6 Extension to the non-Poisson case 



In this section we prove Theorem 2.4. We assume here that all the point processes 
Xn are coupled as described at (3.3), in terms of a sequence ((Xi,Ti), (X2,T2), . . .) 
of independent random elements of M'' x with common distribution k x /x^k. 
Given / G B{R'^) and A > 0, for each m G N we define 



where we set 



m+l,A 



l^m,\) = Ym+l,\ + A 



i.m.X-) 



(6.1) 



1=1 



^m+l,A '■— {f, ^x0^m+l,Tm-^-i; X„i)); 

In this section we shall use the standard notation \\X\\p for the i7-norm E [IXI^]-*^/^ 
of a random variable X, where p > 1. 

Lemma 6.1 Suppose i?(x, T; 'Hk(x)) < oo for n-almost all x G M.'^, and Assumption 
Al or A3 holds, along with Assumption A6' or AT. Let (/i(A))a>i satisfy h{\) / X — > 
as A —> oo. Then with S{x, A) defined at (2.17), 



lim sup 

A^oo A-/i(A)<^<m<A+ft(A) 



EFe,^Frr,,x- (^^«(x)/(x)5(x,«(x))d? 



= 0. 



Proof. We suppress mention of marks in the proof, writing simply X for (X,T). 
Suppose (•^(A))a>i and (m(A))A>i satisfy A - h{X) < e{X) < m(A) < A + h{X). We 
shall show a sequential version of (6.2). To ease notation, we write £ for -^(A), and 
m for m(A), and for Yjn,x and Aj_„ for A^^^^a- 



34 



By using (6.1), expanding and taking expectations, we obtain 



m 

i=i / \ 3=1 j=e+2 



= E [Ye+,Ym+i] + £E[Ai,eYm+i] + iE[Y,+iAi,m] + i{i - 1)E [Ai,,A2,™] 

+{m-e- 1)E A,+2,m] +i{m-i- 1)E [Ai,,A ]■ (6.2) 

We shall establish the hmiting behaviour of each term of (6.2) in turn. First we 
have 

E[Fm^m+i]= / «:(x)(dx) / «:(y)(dy)E[(/,eA(x;A',))(/,eA(y;A'^_i))]. 



Here and below, all domains of integration, when not specified, are IR'^. By Lemma 
3.8, for almost all (x, y) G Fq x Fq, we have 



(/,eA(x; A',))(/,eA(y; A'^_,)) ^ /(x)/(y)e^,(7i.(x),M'^)e^(7i:«(y),M'^). (6.3) 

By (2.16), the variables (/, ^a(x; A'^))(/, ^^(y; '^m-i)) uniformly integrable so 
(6.3) extends to convergence of expectations, and also the limit is bounded, so that 
setting 

71 := ^/(x)E[e^(7i,(.),M'^)]«:(x)(rfx), (6.4) 
we have as A — > oo that 

E[r,+i>;n+i] ^7i- (6.5) 

Next, observe that 

m [Ai,,y;„+i] = m [(/, a(Xi; x,+,) - a(Xi; x,)){f, a(x„+i; x^))] 

J /^(x)cix J K{y)dy J K{w)dwE[{{f,U^;Xr_,) - U^;Xe_^))) 

x(/,6(y;'^m-2U{x,w}))]. (6.6) 

Making the change of variables z = A^/'^(w — x) we obtain 

eE [Ai,,r„+i] = jj «:(x)(cix) J K{y){dy) I «(x + A-^/'^z)(iz 
xE[(/,eA(x;A',lV"''')-eA(x;A',_i)) x (/, eA(y; ^"^-2 U {x, x + A-^/^))]. (6.7) 
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Suppose that Al or A3 holds. By Lemma 3.8, for almost all (x, y) G Fq x Fq, we 
have that 



(/, U^; ^) - U^- X,.,)){f, Ur, U {x, X + A-V'^z})) 



^ /(x)/(y)(e^(K(x),R') - ^Un.i.).^'ML{n.iy).^% (e.s) 



By (2.16) and the Cauchy-Schwarz inequality, the variables in the left side of (6.8) 
are uniformly integrablc. Therefore we have convergence of expectations, so that 
the integrand in (6.7) tends to 

«:(x)^«:(y)/(x)/(y) 
xE [iU^l^.). - ^^(^«(x), [^Un.iy)X)\. (6.9) 

Also, {f\ ^xiy.] Xf^^ ^ ^) — ^a(x; A"^-!)) is zero imless |z| < /2a/-i(x), defined at 
(2.13). Also, by A6' or A7', we assume for some p > 2 and q > 2dp/{p — 2) that 
the moments condition (2.16) holds and we have binomial power law stabilization 
of order q (in A7', since p > 3 the condition q > d{150 + 6/p) ensures that q > 
2dp/{p — 2)). Therefore the Holder and Minkowski inequalities yield 

|E [(/, eA(x; X^:^,'-''"^) - U^; X,.,)){f, Ut, Xra^2 U {x, X + A-^Z-^z}))] | 

< (||(/,a(x;A'r/"''^))||p+ !l(/,a(x; A',_i))||p) 
x||(/,6(y;A'„-2U{x,x + A-i/'^z})||pP[i?,,,_i(x) > |z|]i-2/P 

<const. X (|z|^('-^)/^'Al). (6.10) 

Since q > dp/{p — 2), this is integrable in z. Set 

72:= / «(x)^dx/(x) / dzE[^^{ni^^),R')-^^(n.^^),R% (6.11) 

By (6.7) and the dominated convergence theorem we obtain 

IE [A,,eYm+i] ^ 72 ^ «:(y)rfy/(y)E [e^(7i«(y), M'^)] = 7172. (6.12) 

Next, writing Xi as x, X^+i as y, and X^+i as x + X'^^'^z, we have 

m [y^+iAi,^] ^jj «:(x)dx J K{y)dy j «:(x + X-'/''z)dz 
xE [(/, Ur, '^t-i)) if, eA(x; U {x + A-^Z-^z}) - ^(x; (6.13) 
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By Lemma 3.8 (note that we do not assume that i < m in that result), for almost 
all (x, y) e To X To, we have 

(/, Ur, XF-i)){f, U^; ^1-2 u {x + \-"'^}) - G(x; xi_^)) 

/(x)/(y)e^(H.(y),M'^)(e^(K(,),M'^) - e^(7i:«(x),M'^)). (6.14) 



V 



Using (2.16), we obtain convergence of expectations corresponding to (6.14). Hence, 
the integrand in (6.13) converges to the expression given at (6.9). By a similar 
argument to the one used to establish (6.10), the absolute value of this integrand is 
bounded by a constant times |z|'?(^~p)/pa1, and this is integrable since q > dp/ {p—2). 
Hence, the dominated convergence theorem gives us 

ffi[r,+iAi,J^7i72. (6.15) 

Next, by taking Xi = x, X2 = y, X^+i = x + X'^/'^z and X^+i = y + A~^/'^w, 
we have 

E[Ai,,A2,^] = A-^y «;(x)rfx J K{y)dy J «:(x + A-i/'^z)dz J «:(y + A-^/'^w)dw 

xE [(/, eA(x; XI, U {x + X-y^^z}) - Xl,)) 

X (/, 6(y; u {x + A-v<^z, y + A-v<iw}) - 6(y; u {x + X-'^'z}))]. 

(6.16) 

For almost all (x, y) € Tq x Fq, Lemma 3.8 yields 

(/, eA(x; U {x + A-V<^z}) - eA(x; XU) 
X (/, 6(y; u {x + x-'/\ y + x-'/'w}) - 6(y; u {x + X-'^'z})) 



V 



Also, the quantity on the left is uniformly integrable by the assumption that (2.16) 
holds for some p > 2. Hence we have corresponding convergence of expectations, so 
the integrand in (6.16) converges to 

/(x)/(y)/.2(x)/.^(y)E[e^(?f:(,,),M'^)-e^(7^.W,M'^)] 

Also, we have (/, eA(x; A'/.^Ufx+A-V'^z})) = (/, ^(x; Xl,)) unless i?A/-2(x; {y}) > 
|z| and (/,eA(y; A'^-3U{x + A-V<^z,y + A-V'^w})) = (/, eA(y; ^^-3 U {x + A-V<^z})) 
unless i?A,m-3(y; {x, X + A ^/'^z}) > |w|. Hence, Holder's inequality shows that the 
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absolute value of the expectation in (6.16) is at most 

II (/, 6(x; XI, U {x + A-V'^z}) - G(x; 
X II (/, a(y; U {x + X-'/\ y + A-^Z-^w}) - ^(y; ^^^^3 U {x + A-V'^z})) ||, 

x(P[i?,,,_2(x; {y}) > |z|])(i/2)-V^(P[i?,,^_3(y;{x,x + A-i/'^z}) > |w|])(i/2)-Vf . 

By the assumption (A6' or A7') that moments condition (2.16) holds for some p > 2, 
and that ^ is binomially power law stabilizing of order q > 2dp/{p — 2), this is 
bounded by a constant times 

(|Z|''(2-P)/(2P) Al)(|w|''(2-P)/(2P) Al) 

which is integrable in (z, w). Therefore the dominated convergence theorem applied 
to (6.16) shows that 

£{£-l)E[Ai,eA2,m]^ll (6-17) 
Next, take Xi = x, X^+i = y, X^+i = x + A^^/'^z to obtain 

xE [(/, eA(x; A7_J - eA(x; 6(x; U {x + X-'^'z}) - U^; 

By Lemma 3.8 the quantity inside the expectation tends to zero in probability for 
almost all x, y and all z. Hence its expectation tends to zero as well, since it 
is uniformly integrable by (2.16). Also, the absolute value of this expectation is 
boimdcd by a constant times |z|'^(2-p)/p a 1, by a similar argument to (6.10). Hence, 
dominated convergence yields 

£E[Ai,,Ai,J^O. (6.18) 

Next we have 

E[F,+iA,+i,^] = J K{^)E[{f,U^;X,)){f,U^;X^)-U^;Xm-i))]. (6.19) 

By Lemma 3.8, for almost every x e Fq, we have 

{f,U^;X,)){f,U^;Xj) ^ /(x)2(e^(7^,(,),M^))2 

(/,eA(x; A',))(/,eA(x; A-^-i)) ^ /(x)2(e^(H.(x), M'^))^ 

and using (2.16), wc have the corresponding convergence of expectations so that the 
integrand in (6.19) tends to zero. Also by (2.16), this integrand is bounded, and 
thus 

E [r^+iA^+i,^] ^ 0. (6.20) 
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Next, setting Xi = y, X^+i = x, and X^+i = y + A ^l^i.^ we find that 

xE [(/, Ur, ^^-i) - 6(y; Xi-i)){f, U^; u {x + A-v<iz}) 

(6.21) 

By Lemma 3.8, for almost all (x, y) G Fq x Fq and all z, as A — > oo we have 

(/, Uy\ xU) - Ur, Xe-i)){f, Ca(x; xy_, u {x + x-'/'z}) - eA(x; 0, 

so that the quantity inside the expectation in (6.21) tends to zero in probability 
and by (2.16) it is uniformly integrablc. Hence the integrand in (6.21) tends to 
zero. Also, by a similar argument to (6.10), the absolute value of this integrand is 
bounded by a constant times |z|^(^~*')/*' A 1, which is integrable since q > pd/{p — 2). 
Thus, the integrand in (6.21) is bounded by an integrable function of (x, y, z) so the 
dominated convergence theorem shows that 

eE [Ai,,A,+i, J ^ 0. (6.22) 

Next, write X^+2 as x, X^+i as y, and X^+i as x + A'^/'^z, to obtain 

E [Yi+^Ai+2,m] = A-' J «:(x)dx J K{y)dy j «;(x + A-^/''z)dz 

xE [(/, Ut, ^d) if, 6(x; U {x + X-'/'^z}) - 6(x; (6.23) 

By a similar argument to (6.10), the absolute value of the expectation inside the 
integral is bounded by a constant times |z|^(^~*'^/*' A 1, which is integrable since 
q > dp/{p — 2). Therefore, the triple integral in (6.23) is bounded, and since m — 
£ — 1 — o(A), it follows that as A — > oo we have 

{m-e- 1)E [F,+iA,+2,m] ^ 0. (6.24) 

Next, put Xi — X, X^+2 = y, X^+i = X + A'^^'^z, and X^+i = y + A~^/'^w, to 
obtain 

[Ai,£A^+2,m] = A~^ J K(x)dx J K{y)dy J dz J dw 

x«:(x + A-i/'^z)K(y + A-i/^w)E [(/, ^(x; A't'l'"''') " U^'^ ^^-0) 

X (/, 6(y; ^"^-3 u {x + x-'/'z, y + x-'/'w}) - Ur, x^-, u {x + x-'/'^z}))]. 

(6.25) 



E 
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By the argument used in dealing with E [Ai^^A2,m] above, the absolute value of the 
integrand in (6.25) is bounded by a constant time's (|z|'?(2-p)/(2p) a1)(|w|'?(2-p)/(2p) a1), 
and hence the integral in (6.25) is bounded. Since £(m — £—!) = o(A^), this shows 
that 

e{m-e- 1)E [Ai,^A^+2,m] ^ 0. (6.26) 

We have obtained limiting expressions for the nine terms in the right hand side of 
(6.2), namely (6.5), (6.12), (6.15), (6.17), (6.18), (6.20), (6.22), (6.24) and (6.26). 
Combining these, we find that as A ^ oo, E [F^Fj„] converges to (71 + 72)^- Since 
the choice of £(A),m(A) is arbitrary, we then have (6.2). □ 

Lemma 6.2 Suppose the assumptions of Theorem 2.4 hold. Suppose h{X) is defined 
for A > 1 satisfying h{X) > and h{X)/X — > as A — > 00. Then F^^x defined at 
(6.1) satisfy 

hm sup (E[F^J)<oo. (6.27) 

A->-oo \-h{\)<m<\+h(\) 

Proof. We abbreviate notation as in the preceding proof. Note that our assumptions 
(in particular A6' or A7') imply that for some p > 2 and q > 2dp/{p — 2), (2.16) 
holds and ^ is binomially power-law stabilizing of order q. 

Let m = m(A), A > 1, be defined to satisfy m(A) ~ A as A — > 00. By a similar 
expansion to (6.2), we obtain 

E [F^J = E [y^+i] + 2(m - 1)E [y„+iAi,„] + m(m - 1)E [Ai,^A2,„] + mE [AjJ. 

We consider these terms one by one. First, E is bounded by (2.16). Second, 

setting Xi = X and X^+i = x + X'^^'^z we have that 

2mE Ai,^] = — / rix / dzK(x)K(x + A~^/'^z) 
xE [(/, eA(x + X-'/'z; A-^.J) (/, eA(x; X^t^^) - U^; X^-i))] . 

Since (/, 6(x; A'^lt '^'^)) = (/, Ca(x; A'^-i)) when i?A,m-i(x) > |z|, use of Holder's 
inequality, followed by (2.16) and the binomial power-law stabilizion, shows that the 
absolute value of the expectation in the integrand is bounded by 

II (/, eA(x + x-'/\ x-_,))u (/, eA(x; x-t^'"! - a(x; x^^^m, 

x(P[i^A,„^-l(x) > |z|])l-(2/^') 

<const. X (jzl^^^-fVP Al), (6.28) 
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which is an integrable function of z. This shows that 2mE [Fm+iA^] is bounded. 
Next, take Xi = x, X^+i = x + A~^/'^z, and X2 = x + A~^/'^(z + w), to obtain 

E[Ai,^A2,r„] =^~^y" /^Wc^x J «;(x + A-^/'^x)dz J k{x + X'^^'^iz + w))dw 

xE [(/, U^; X-tV"^^^""^ U {X + A-V<^z}) - U^-, KtV'^'^^^)) 
X (/, + A-V<^(z + w); X-_, U {x + A-V<^z}) - + A-V<^(z + w); _,))]. 

(6.29) 

Inside the expectation, the first factor is zero if -RA,m-2(x; {x + A^^/'^(z + w)}) < |z|, 
and the second factor is zero if i?A,m-2(x + A~^/'^(z + w); {x}) < |w|. Hence by 
Holder's inequahty, the assumption (2.16), and the assumption of binomial power- 
law stabilization of order g, the expectation inside the right side of (6.29) is bounded 

by 

II (/, eA(x; X-tV^'^^-"-^ U {X + A-V<^z}) - U^- X-tV'^^^^^Mv 
X II (/, + A-V<^(z + w); X^_, U {x + A-V<iz}) - a(x + + w); Xl_Mv 

x(P[i?A,m-2(x; {x + A-^/'^z}) > |z|])(i/2)-Vp 

x(P[i?A,m-2(x+A-V'^(z + w);{x}) > |w|])(V2)-l/p 

< const. X (|z|«(2-P)/(2p) /\ 1) X (|w|'?('-^')/(2^') A 1), 

and since q > 2dp/ (p — 2), this uniform bound is integrable in z, w. This shows that 
m(m — 1)E [Ai^toA2,to] remains bounded. 

Finally, take Xi = x and X^+i = x + A~-'^/'^z, to obtain 

mE[AU = jjdxj cizK(x)K(x + A-'/'^z) 

xE[(/,{A(x;A'^ir^''^) -{a(x; A'„_i))']- (6-30) 

Since the quantity inside the expectation is zero unless i?A,m-i(x) > |z|. Holder's 
inequality shows that this expectation is bounded by 

(E[|(/,6(x;A'^tr^^'^^) -6(x; A'^_,))n)VP(P[i?,,_,(x) > |z|])^-^/^ 

< const. X {\z\i^^-PyP Al), 

which is integrable in z. Hence, mE [A^ ^^j is also bounded. □ 

Proof of Theorem 2.4- Suppose ||k||oo < 00 and k has bounded support. Suppose 

^ is k(x)— homogeneously stabilizing at x for K,-almost all x G M.'^, satisfies Asuump- 
tion A4 or A5, and also satisfies A6' or A7'. Let / e B{M.'^), and assume either Al or 
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A3 holds. Suppose (A(n))n>i is a (0, oo)-valued sequence with |A(n) —n\— 
as n ^ oo. 

Let Hn := (/, ^i(n) n) •= (/) f^x{n))- "^^^ ^^^^ proof, assuHie that for all n, 

Xn is given by (3.3) and that V\i^n) is coupled to Xn by setting V\{n) — U^]^{(Xj, Tj)}, 
with an independent Poisson variable with mean \{n). Let 



a :- 



J /(x)5(x, fi;(x))fi;(x)(ix. 

First we show that as n — > oo, 

E [{n-'/\H'^ -H^- (iV, - n)a)f] ^ 0. 
To prove this, note that the expectation in the left hand side is equal to 

n"^ ^ E [((/, vx{n),m - ^\{n),n) - (m - n)af] P[Nn = m] 

m:\m- X{n)\<n^/* 

+n-^E \{H^ -Hn- {Nn - n)af l{|iV, - A(n)| > n^/^} 



(6.31) 



(6.32) 



Let £ > 0. By (6.1) and Lemmas 6.1 and 6.2, there exists c > such that for large 
enough n and all m with A(n) < m < A(n) + n^/^. 



E [((/, mn),m-'^\{n),n)-{m-n)af 



E 



m—l 



e=n 



< £(m— n)^+c(m— n). 



where the bound comes from expanding out the double sum arising from the expec- 
tation of the squared sum. A similar argument applies when A(n) — n^^^ < m < n, 
and hence the first term in (6.32) is bounded by the expression 

n-^E [e{Nn - nf + c\Nn - n\] 
< n-\e{X{n) - nf + eE {{N^ - K"^)f\ + cE [|Ar„ - A(n)|] + c|A(n) - n|) 
< n-^{e{\{n) - nf + e\{n) + cA(n)^/2 ^ ^|^(^) _ ^|^^ 

and so, since e is arbitrary, the first term in (6.32) tends to zero. 

Our assumptions (A6' or A7') include the moment bounds (2.2) and (2.16) for 
some p > 2. Hence, choosing p' G (2,p) we can apply Lemma 4.3 of [19], taking 
p\ = A^/^^*^-* in that result, to bound the L^' norm of the contribution to H' from 
points in a cube of side \(n)~^^^'^'^^ by 0(A(n)^^^^''^^^^). The number of such cubes 
intersecting supp(/t) is 0(A(n)^/*^^'^-'), so that by Minkowski's inequality we obtain 



\Hn\\p' 



(6.33) 
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Also, the value of Hn is the sum of contributions from the n points of Xn, and 
by (2.16), the norms of each contribution are bounded, so that by Minkowski's 
inequality ||i^n||p = 0{n) so that = 0{n) also. Moreover, ||A^n||p' — 0{n). 

Combining these facts with (6.33), we may deduce that 

\\K -Hn- {Nn - n)a\\,, = 0{n^'^+'^/^'^^). 

Hence, by Holder's inequahty the second term in (6.32) is bounded by a constant 
times n'~-^n^^^+^)/^(P[|A^n-A(n)| > n^/^])^~(^/^'\ which tends to zero (see e.g. Lemma 
1.4 of [13]). This completes the proof of (6.31). 

Set (7^ := Jp /(x)^V^^(x, «;(x))«;(x)(ix. and set t| := cr^ — a^. By Theorems 

2.2 and 2.3 we have as n ^ cx) that Var(ij;j ^ a"^ and n-^''^{H'^ - Eif^) 
A/'(0,(T^). Using (6.31) and following page 1620 of [16] verbatim, we may deduce 
that lim„^oo n~^Var(if„) — > Tj and also 

^"'/'a^i(n),n) =^"'/'(^n-E//„) ^AA(0,r|). (6.34) 

Since t| is the limiting variance in (2.18), we thus have (2.18). 

Suppose that /i, . . . , /fe are in BiW'') or that /i, . . . , are in B{W^) and Al holds. 
If ai, . . . , afc are real constants, by (6.34) we obtain convergence of Yl!i=i ^~^(^i{fii ^\{n),n) 
in distribution to the centred normal with variance 

^^aittji / fi{x)fj{x)V^{x,K{x))K{x)dx 
i=i 3=1 

/i(x)5(x, fi;(x))fi;(x)(ix / /j(x)(5(x, fi;(x))fi;(x)(ix 



with V^^(x, a) given by (2.15). The convergence in distribution follows by the 
Cramer- Wold device. □ 



7 Applications 

Many examples and applications of the general theory are described in [3, 16, 17, 
18, 19], and here we discuss only a few. The examples we consider have translation- 
invariant ^. There are interesting potential applications of the theory with non 
translation-invariant ^ to topics in multivariate statistics such as nonparametric 
density estimation and nonparametric regression [5, 4, 7], but these are not easy 
to describe briefly in an already lengthy paper. The first example discussed here 
illustrates the application of the Law of Large Numbers in Theorem 2.1. 
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7.1 Voronoi coverage 

For finite A' C M'' and x e A", let C(x; A") denote the Voronoi cell with nucleus x for 
the Voronoi tessellation induced by X. Given a density function k on R'^, and A > 0, 
let Vx denote a Poisson point process in with intensity function \k{-). With a 
view to potential applications in nonparametric statistics and in image analysis, 
Khmaladze and Toronjadze [10] ask whether, for an arbitrary bounded Borel set 
A C R*^, the total volume of bounded cells C(x; Vx) with nuclei at points x e P;^n A 
converges almost surely to the Lebesgue measure of as A — > oo. They also ask 
whether 

\AA U^^r,nAC{x;Vx)\^0, 

where | • | here denotes Lebesgue measure. They answer these questions affirmatively 
for the case d = 1 only. If one is satisfied with convergence, one can use our law 
of large numbers (Theorem 2.1) to answer the first question affirmatively for general 
d, and to partially answer the second question also. We also have corresponding 
results for X^- It does not seem possible to achieve these results using only the 
results of [18] or [3]. 

To put these questions in our framework, define the set V{'x; X) to be the set 
C(x; X) if this set is bounded, and to be the empty set otherwise. For finite A C M"* 
and X e A", let ^(x; A, ■) be the restriction of Lebesgue measure to F(x; A"), and 
let ^*(x; X, ■) be the corresponding point measure, defined at (2.4) Note that in this 
case, ^ is translation- invariant and points do not carry marks. 

Our choice of ^ has the homogeneity property of order rf, which says that 
^(ax;aA,M'^) = a'^{(x; A, M'^) for any a > 0. Hence, for Borel A C R"^, using 
(2.1) we have 

A-Vf(A)=A-i 6(x;Pa,M'^) = A-i e(A^/'x;AV'^P,,R<^) 

= E e(^;^A,M'^)= Yl inx;PA)|. 

Similarly, 

xeAnAVi 

By arguments in [16], the measure ^ satisfies R{0;'Hx) < oo almost surely for all 
A > 0. Also, 

m^inx,^')]^i/x 

since the average volume of Voronoi cells in a homogeneous Poisson process of in- 
tensity A must be 1/A (one can use Theorem 2.1 to show this rigorously). Moreover, 
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by arguments in [16] the measure ^ satisfies the moments conditions (2.2) and (2.3), 
for example if k is supported by the unit cube in and bounded away from zero 
on its support. Hence, under these conditions on k, by setting / to be the indicator 
function of 1^, we can apply Theorem 2.1 to deduce that 

A-Vr(^)-^ / (l/«;(x))«;(x)dx= |^nsupp(K)| 

J j4nsupp(/t) 

and likewise if \{n) n as n ^ oo, then n uj^^^^ ^{A) — ^ \A fl supp(k)| for 
arbitrary Borel ^4 C R*^. This answers the first question raised above, in the 
sense. 

If C(x, X) is bounded, then with / = 1^ we have 

(/, U^; X) - e:(x; X)) = A(|C(x, X) n A\ - |C(x, A')|l^(x)) 
= A(|C(x, X) n A|1ac(x) - |C(x, X) \ ^|1a(x)) 

so that if there are no unbounded Voronoi cells having non-empty intersection with 
A, then 

\AA U^eP.nA C(x; Vx)\ - A-^ Yl ^^(^; ^a) - r(x; Pa)) I 

which converges in to zero by Theorem 2.1. In cases where the closure of A is 
contained in the interior of supp(/t), the probability of there being any unbounded 
Voronoi cells intersecting A tends to zero, so the above partially answers the second 
question raised above in an sense. However, we do not fully deal here with sets 
touching the boundary of supp(«;). 

7.2 Germ-grain models 

Germ-grain models are a fundamental model of random sets in stochastic geometry; 
see for example [8, 12, 23]. In the germ-grain model, a random subset of is 
generated as the union of sets (X^ + Ti) where {Xi} (the germs) are the points of a 
point process, and {Tj} (the grains) are independent identically distributed random 
compact subsets of W''. Our results can be applied to obtain hmit theorems for 
random measures associated with germ-grain models, in the case where the point 
process of germs is V\ or X^, and where the grains are scaled by a factor of A"^/'' 
as A — > cxD. 

Let Ai denote the space of compact subsets of M*^. For finite X G M.''- x A4, and 
A > 0, set X^ := {{x,X-^/H) : (x,t) e X}, and set 

SA(A') = U(x,t)eAr(x + A-n). 
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When X is Vx or Xn, the set Sa(A') is a germ-grain model with germs given by a 
Poisson process or binomial process and grains scaled by a factor of A~^/'^. We can 
apply our general results to the volume measure of Sa(Pa) (i-e., the restriction of 
Lebesgue measure to Sa(Pa)) and the surface measure of Sa(Pa) (i-e., the restriction 
of {d — l)-dimensional Hausdorff measure to the boundary of Ex{Vx)), and hkewise 
for Xn- 

For t E M., i.e. for t a compact set in M*^, let \t\ :— max{|z| : z E t} and let \\t\\ 
denote the Lebesgue measure of t. 

Theorem 7.1 Suppose that for q = 1 or q — 2, and for some p > q, we have 
E[\\T\\P] < oo. Then for f G B{M.'^) the integral J„^^^^^ /(x)(ix converges in L'^ to 

a finite non-random limit. If A(n) ~ n as n — > oo, the integral /-.^^ ^^.^^ ^/(x)cix 
converges in to the same limit. 

We sketch the proof. For finite X G M.'^ x Ai , let 7r{X) denote the projection of 
X onto Mf^, i.e. the subset of Mf^ obtained if we ignore the marks carried by points 
of X. Also, for each x G 7r(A') let T(x) denote the mark carried by x, i.e. the value 
of t such that (x, t) G X. For y G Si (A") = U(x,t)GA'(x + t), let Nx{y) denote the 
nearest point x G 7r(A') to y such that y G x + T(x) (in the event of a tie when 
seeking the 'nearest point', use the lexicographic ordering as a tie- breaker) . Take 
^{^,t]X,-) to be the restriction of Lebesgue measure to the set of y G x -|- i such 
that X = Ux{y)- 

Then, since Nxij) is unique for each y G Si(A'), 'Yli{:x.t)ex ^(^' ') precisely 
the volume measure of Si(A'). Also, { is translation-invariant, so by (2.1), 

a(x, t- X, A) = e(A'/'^x, t; X^^'^X, X'/'^A) = X^{x, X'^'% X^, A), 

so that Yli^^ t)ex Ca(x, t; X, ■) is the volume measure of Si (A''*'), which is the same 
as the volume measure of Ex{X). Hence, with this choice of ^, we have that 

A-Vi(c?x) = ls,(P,)(x)dx; (7.1) 
>^M,nid^) = ls.(Af„)(x)dx. (7.2) 

The measure ^(x, t; A", •) is supported by x + t, and this measure is unaffected by 
changes to X outside B2\t\{'X-) x M.. This is because for any y G x -|- i it is the case 
that |y — AGf(y)| < 1^1 so by the triangle inequality, Mx{y) lies in _B2|t|(x). Hence, 
2\t\ serves as a radius of stabilization, which is almost surely finite since M. is the 
space of compact subsets of W^. Also, .^(x, t; A, R'^) is bounded by and the 
conditions (2.2) and (2.3) follow. 

We can then apply Theorem 2.1 to obtain the result. The measure ^(x, t; X, ■) 
satisfies Assumption A2, so we can take test functions in B{W^). The limit is given 
by the right hand side of (2.5). 
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By applying Theorems 2.3 and 2.4 to the above choice of ^, we obtain the fol- 
lowing result on convergence to a Gaussian field for the volume measure on SxiVx) 
or on SA(n)('^n)- We state the result in terms of the measures and z/^^, which 
translate into statements about the volume measures by (7.1) and (7.2). 

Theorem 7.2 Suppose k is bounded and has bounded support. Suppose for some 
p > 3 that E[||T||^] < oo and there exists C > and q > d{150 + 6/p) such that 
P[\T\ > s]< Cs-1 for all s. 

Then with ^ as given above, the finite dimensional distributions of the random 
field X~^^'^{f,]j^^) , f e S(M'^), converge to those of a centred Gaussian random fi,eld 
with covariances given by Jj^^ /i(x)/2(x)y^(x, /t(x))(ix, with given by (2.15) and 
^ as defined above. Likewise, if |A(n) — n\ = 0{n~^^'^), then the finite dimensional 
distributions of the random field „), / G BiW^), converge to those of 

a centred Gaussian random field with covariances given by the right hand side of 
(2.19). 

Theorem 7.2 adds to the results for germ-grain models in ([3], Section 3.3) in 
several ways. In particular, in [3] it is assumed that the distribution of \T\ is sup- 
ported by a compact interval, whereas here we need only power-law decay of the tail 
of this distribution. Also, in [3] the term 'volume measure' is used in a non-standard 
way to refer to an atomic measure supported by the points of X. Our usage of the 
terminology 'volume measure' seems more natural, and is also in agreement with 
the standard usage found, for example, in [9, 23]. It is not clear whether the general 
results in [9] can be applied to the volume measure. 

We now consider the surface measure of SxlX). We assume here that with 
probability 1, each grain is a finite union of bounded convex sets. For (x, t) G X, 
let (x -I- ty and 5(x -|- t) denote the interior and boundary, respectively, of the set 
:x. + t. For z G U(x,t)GAr(x -|- t)°, let J^xi'^) denote the closest point x G to z 

such that z G (x + T(x))°, using lexicographic ordering as a tie-breaker. 

Define the set NC{x,t) (i.e., the points for which x is the 'nearest covering' 
germ) by 

NC{x, t) = {z G (x + t)" : x = Ar^(z)} 

which is the set of points interior to x -|- i which are not covered by a set with germ 
closer than x. Define the set CC(x, i) (the points of d{'x. + t) with 'closer cover') by 

CC(x, t) = U(y,„)e;,\{(x,i)}{z G 5(x + t) n (y + m)° : y = A6t(z)}. 

Let us take ^(x, t; X, •) to be the following signed measure: 

• Let ^+(x, t; X, ■) be the restriction of {d — l)-dimensional Hausdorff measure 
to d{x + t)\CC{x, t). 
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• Let ^ {x.,t;X, •) be the restriction of {d — l)-dimensional Hausdorff measure 
to the set 

NC{x, t) n (U(,,„)e;t\{(x,t)}5(z + u) \ CC{z, u)) . 

• Let ^(x, t; X, ■) be the signed measure ^''"(x, t] X, ■) — ^~(x, t] X, ■). 

The signed measure ^(x, t; X , ■) is supported by x+t and is unaffected by changes 
to X outside i?2|t|(x) x M.. Achieving this is the purpose of the definition of ^ used 
here, since it ensures \T\ serves as a radius of stabihzation. 

We assert that J^^^ ^^g;^, ^(x, t; X , •) is precisely the surface measure of Si (A"). To 
see this, suppose z hes on the surface of x + t, but is covered by some other {y + u)° 
with (y,M) G X (take the closest such y to z). If |y — z| < |x — z|, then z G CC(x,t) 
so that ^"""(x, t] X, dz) = ^~(x, t; X, dz) — 0. If |y — z| > |x — z|, then z ^ CC(x, t) 
so that z e NC{y,u) n 9(x + t) \ CC{x.,t), so that ^"""(x, A", dz) is the surface 
measure of 5(x + 1), and ^~(y, u; X, dz) is also the surface measure of 9(x + 1), and 
these cancel out. If also z G (w + v)" with {w,v) G A" and |w — z| > |y — z|, then 
z ^ ^^(^(w,^) so that ^~{w,v;X,dz) = 0. 

Since ^ is translation-invariant, by (2.1) we have 

eA(x, t; X, A) = ^{X^/'^x, t; X^/'^X, X^Z-^A) = X^'^-^^/'^^{x, X-^'H; X^, A) 

and hence, X^^~'^^/'^n\ is the surface measure of 'Ex{V\)^ while X^^'''^^/'^v\^ is the 
surface measure of Sa(A'„). 

To apply our general results here, we need the moments conditions such as (2.2) 
to apply to both the positive and negative parts of the measure ^. We write \dT\ 
for the [d — l)-dimensional Hausdorff measure of the boundary of T. Clearly this 
is an upper bound for ^+(x, T; A',R"'). To estimate the negative part, observe that 
all contributions to ^~{yi,T]X) come from the boundaries of sets associated with 
germs within distance at most 2171 from x. Hence, we have that 

where Xi are independent copies of |9T|, and N is Poisson with mean 9{2\T\Y, 
with denoting the volume of the unit ball. By Minkowski's inequality, the above 
expectation is bounded by E [A'"*']E [Ic^Tl**], and hence the condition 

E[|T|*]E[|9rn < oo (7.3) 

suffices to give us all the moments conditions (2.2), (2.3), (2.7) and (2.16), for both 
the positive and the negative parts of ^. 
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Thus, with this choice of ^ we may apply Theorem 2.1 (with test functions 
/ e B(M.'^)) if for g = 1 or 5 = 2 we have (7.3) for some p > q. We may apply 
Theorems 2.3 and 2.4 (again with test functions / G B(M.'^)), either if (7.3) holds 
for some p > 2 and P[\T\ > r] < Ce~'^^'" for some C > and all r > 0, or if for 
some p > 3 and some q > d{150 + 6/p), (7.3) holds and P[\T\ > r] < Cr~'^ for some 
C > and aU r > 0. 

7.3 Random packing measures 

The random packing measures discussed in Section 3.2 of [3] are obtained by particles 
(typically balls) being deposited in space at random times, according to a space-time 
Poisson process. Particles have non-zero volume and (in some versions of the model) 
may grow with time, but deposition and growth are limited by an excluded volume 
effect. In [3], the measures associated with these packing processes are obtained 
as a sum of unit point masses, with one point for each particle. As in the case of 
germ-grain models, it is quite natural instead to consider the the volume measure 
associated with the random set obtained as the union of particles (balls), or even 
the surface measure of this random set. The setup of this paper enables us to do 
this, but we do not give details. 
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